


% 

" -r. 






NAVAL POSTGRADUATE SCHOOL 






AN INVESTIGATION OF SEVERAL BRANCHING 
FUNCTIONS IN A BRANCH AND BOUND ALGORITHM 
FOR THE CHROMATIC NUMBER PROBLEM 

by 

Ronald Ernest Rautenberg 
December, 1980 

Thesis Advisor: Douglas R. Smith 

Approved for public release; distribution 
unlimited 



T197858 



seCUNITV CLASSiriCATION or this r Aoe (W»m> Oma Entmra^ 



^UOir y». 



REPORT DOCUMENTATION PAGE 



1. ne^oAT NUMacn 



2. GOVT Accession NO. 



4. T|TL£ (•fid Subtitlm) 

An Investigation of Several Branching 
Functions in a Branch and Bound Algorithm 
for the chromatic number Problem 



' V READ INSTRUGVriONS 
BEFORE COMPLETlSdS^ORM 



3. nect^l Cn T*S CAT ALOG number 



5. TYRE OF REPORT 4 PERIOD COVERED 

Master's Thesis; Dec 80 



S. PERFORMING ORC. REPORT NUMBER 



7. AUTMORr#> 



i. CONTRACT OR GRANT NUMBERr*) 



Ronald Ernest Rautenberg 



f. PERFORMING ORGANIZATION NAME ANO ADDRESS 

Naval Postgraduate School 
Monterey, California 93940 



10. program ELEMENT. PROJECT TASK 
AREA 4 WORK UNIT NUMBERS 



CONTROLLING OFFICE NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, California 93940 



>2. nePOAT DATE 

December , 1980 



13. NUMBER OF pages 



12± 



U. MONITORING AGENCY N AME * ADDRESS<f# dUimerntt tr«m CamtnUIng 

Naval Postgraduate School 
Monterey, California 93940 



It. SCCUPITY CLASS, (el Ihle ripert) 

Unci ass i fi ed 



ts«. declassification/ DOWN grading 

SCHEDULE 



tS. DISTRIBUTION STATEMENT fol R«p«rO 



Approved for public release; distribution unlimited 



17. DISTRIBUTION STATEMENT (of thm mbotrmct mntotod in Btock 20, If dlffotont from Mopori) 



IS. supplementary NOTES 



If. KEY WORDS (Contin%»o on rowotoo oldo if nocoooonr arW idonUty df Mock nuMS«r> 

chromatic number, graph coloring, k-chromatic, 
branch and bound, vertex ordering, articulation points 



20. abstract (Contiwmao on rowmmo cido it nocoomorj ond idoniitf ky kloek mmikor) 

The chromatic number problem is to determine the minimum number 
of colors to assign to the vertices of a graph such that no connect? 
ed vertices are assigned the same color. This paper presents a 
branch and bound solution to the chromatic number problem and in- 
vestigates five different branching functions. Additionally, a 
method of coloring very sparse graphs is presented which divides a 

graph into biconnected components and reduces the time required to 
color the ora oh. 



DD 



FORM 

I JAN 73 



1473 EDITION OF I NOV SS IS OBSOLETE 
S/N 0103*014-6«01 I 



1 



SECURITY CLASSIFICATION OF THIS PAGE Dofm ktuocod) 



Apopoved for public release? distribution unlimited 



An Investigation of Several Branching Functions 
in a Branch and Bound Algorithm for the 
Chromatic Number Problem 



by 



Ronald E. Rautenberg 

Lieutenant Commander^ United States Navy 
B.S.» University of Washington# 1970 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN COMPUTER SCIENCE 

from the 

NAVAL POSTGRADUATE SCHOOL 
December 1900 



ABSTRACT 



The chromatic number problem is to determine the minimum 
number of colors to assign to the vertices of a graph such 
that no connected vertices are assigned the same color. 
This oaper presents a branch and bound algorithm for the 
solution to the chromatic number problem and investigates 
five different branching functions. Addi t i ona 1 1 y » a method 
of coloring very sparse graphs is presented which divides a 
graph into biconnected components and reduces the time re- 
quired to color the graph. 
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INTRODUCTION 



Informally^ a graph is a collection of points^ called 
nodes or vertices# some of which may be connected by lines# 
called edges. A solution to the graph colorability problem 
specifies a color for each node such that any two nodes 
which are connected by an edge oo not have the same color. 
The optimal colorability problem is to assign colors# as in 
the graph colorability problem# with the restriction that 
the fewest number of colors possible is used. 

For centuries# mathematicians have studied graphs and 
their characteristics with the awareness that many engineer- 
ing and mathematical problems could be modeled by a graph. 
Then questions about the engineering problem could be 
answered by studying the c ha r ac t e r i s t i c s of the graph and 
using knowledge derived from graph theory. 

One of the earliest and best documented graph colorabil- 
ity problems dates back to 1852# when a student of De Morgan 
named Francis Guthrie# proposed that 4 colors was all that 
was needed to color a mao such that no two bordering coun- 
tries are the same colorlRef. 11. Note that this problem is 
transposed into a graph colorability problem by using a node 
for each country# and connecting 2 nodes with an edge if the 
two corresponding countries share a common border. 
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In fact/ beginning in the 19th century/ many of the en- 
gineering problems in areas such as scheduling/ communica- 
tions/ transportation/ electronics/ and chemistry have been 
studied through the use of graph theory. More recently/ the 
social/ biological/ and environmental sciences have relied 
on graph theory to investigate specific problems. In the 
decade of the seventies/ computer science and graph theory 
have been closely intertwined as computer scientists real- 
ized that many of their difficult problems can best be un- 
derstood by graph modeling. One of the most notable is the 
scheduling problem where two or more processes must compete 
for common resources. Ways to efficiently schedule these 
processes are better understood when the system is modeled 
by a graph. 

The intent of this research has been to take a single 
problem from graph theory/ namely optimal graph colorabili- 
ty/ and to investigate the use of the digital computer to 
find optimal colorings for various graphs. Though a seem- 
ingly trivial problem on a computer/ it is complicated by 
the fact that in the worst case/ all known methods of find- 
ing solutions to problems of this type require exponential 
time. That is/ as the size of the graph gets larger/ the 
time required to find the optimal coloring grows exponen- 
tially. 
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Tyoically/ the researcher investigating an exponential 
time problem, whose size is very large, is faced with one of 
two choices: 1) he must accept a solution which is not 
necessarily optimal/ but some approximate or 'close* solu- 
tion, or 2) he must relax some of the restrictions on the 
problem. Most of the research done on the graph colorabili- 
ty problem has been in the area of finding approximation al- 
gorithms, that is, algorithms which attempt to find a rela- 
tively few number of colors without undue computational ef- 
fort being expended tRe f s . 2, 3,^1. The primary thrust of 
this research has been to investigate methods of finding 
only optimal colorings, and to study the effects of dif- 
ferent procedures on the time and soace (computer memory) 
required by each. 

Specifically, the branch and bound technique is used to 
find the optimal coloring of a large number of randomly gen- 
erated graphs, and the branching function is varied to 
determine the relative value of different heuristic branch- 
ing functions. 

In summary, this research encompasses three related, yet 
diversified disciplines. From mathematics comes graph 
theory; from ooerations research comes the branch and bound 
techniques; and from computer science comes problem complex- 
ity, Background information is provided for these three to- 
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pics in the first part of this report. Following that is a 
discussion of the research done on the computer/ including 
the algorithms used/ the branching functions investigated/ 
and the graph generation techniques. Finally/ there is a 
discussion of the results obtained/ the difficulties encoun- 
tered/ and the ideas and notions derived during the course 
of the research. 
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II 



BACKGROUND 



A. graph theory 

While not attemoting in any way to fully cover the field 
of graph theoryf an introduction to the notions and termi- 
nology which relate to the coloradility problem is necessary 
for an understanding of the discussion which follows. The 
definitions used here are standard in most of the literature 
on graph theory^ so the reader with some background in the 
field can skip this section without any loss of continuity. 

Formally/ a graph G = (V/E) consists of a finite/ 
nonempty set of nodes or vertices V/ and a finite/ possibly 
empty set of vertex pairs called edges E. For the purooses 
of this research the set of edges is restricted to not con- 
tain any elements of the form (x/x). That iS/ no edge can 
exist from a vertex to itself. In the case where the edges 
have direction/ they are defined by an ordered pair of ver- 
tices and the graph is called a DIGRAPH. This research has 
been strictly limited to non-directed graphs. 

The general graph colorability problem is formalized as 
follows: Given graph G and integer k/ does there exist a 
function R: V -> I where I = (1/2/. ../k> such that for all 
(x/y) S E/ R(x) * R(y)? The optimal colorability problem is 
to find the smallest integer k such that a function R ex- 
ists. That smallest integer found is called the chromatic 



number of the graph G. If R(x) = R(v) then vertices x and y 
are said to belong to the same color class. Note that the 
function R partitions the vertices of the graph into k color 
classes. 

Some additional terms relating to graphs are: 

Adjacent: Vertex x is adjacent to vertex y if an edge ex- 

ists between them# that is (x>y) 6 E. 

Adjacency matrix: A two dimensional/ n x n matrix# used to 

describe any graph. The indices into the matrix are 
the vertices of the graph# and the elements are a 1 
if an edge exists between the two vertices# and a 0 
otherwise. For the graphs used in this research the 
diagonal elements are all zeroes and the matrix is 
reflective about the diagonal. The adjacency matrix 
is a very common# and easy to implement method of 
representing a graph in a computer. 

Degree: The degree of a vertex is the number of edges in- 

cident upon it# or the number of vertices which are 
adjacent to it. The degree of vertex x is easily 
calculated from the adjacency matrix by simply sum- 
ming the I's in row x or column x. If the elements 
of the adjacency matrix are denoted a(i#j) and the 
degree of vertex x is denoted d(x) then 

d(x) = XaCx^i) 
i 
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Path: A oath is said to exist fpoin vertex x to another ver- 

tex y if one can start at vertex x and follow the 
edges of the graph and reach vertex y. Formally^ a 
path exists if there exists an ordering of some of 
the vertices of G {vjfV 2 /»../v^> such that vj = x f 
v^ = y and for 1 <= i <= n-W (v^/v^-^j) 6 E. 

Connected Graoh: A graph is connected if at least one path 

exists from each vertex to every other vertex in V. 

Doubly connected graoh: A graoh is doubly connected if for 

every vertex x and every vertex y at least 2 oaths 
exist from x to y and no other vertices are common 
to both oaths. This is also called a biconnected 
graoh . 

Articulation point: A vertex v is an articulation point of 

graph G if there exists some vertex x and some ver- 
tex y such that every oath from x to y includes ver- 
tex V. Note that a doubly connected graoh has no ar- 
ticulation ooints. 

Complete graph: A graph G is complete if all possible edges 

are present. That is G is complete if for all ver- 
tices X and y» (x/y) 6 E. 

Cligue: A digue is a subset of a graph such that every 

vertex in the clique is adjacent to every other ver- 
tex in the cli que . 
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0. BRANCH AND BOUND 



Many oroblems which deal with searching for a solution 
or set of solutions satisfying some constraints can be 
solved using the branch and bound technique. It is esoe- 
cially useful in solving minimization (or maximization) 
oroblems and has been successfully used in constrained op- 
timization oroblems since the late 1950's. Although there 
are several ways of describing branch and bounds the basic 
idea is to split the solution space (branch) and place a 
limit# or lower bound# on the optimal cost of the problem 
limited to each subset of the solution space. Those subsets 
whose optimal costs do not exceed the cost of some known 
solution are then repeatedly divided and bounded until a 
solution is found whose actual cost is no greater than the 
lower bound on any of the subsets. This solution is thus 
the optimal one. The power of branch and bound comes from 
its ability to leave unexplored those supsets which are 
known not to contain the optimal solution. 

In order to use the branch and bound technique the prob- 
lem must have an expressable set of solutions and a cost 
function COST() which maps the solutions into nonnegative 
integers. Usually# the problem calls for finding that solu- 
tion s for which COST(s) is a minimum. Sometimes all solu- 
tions may be desired whose cost is no greater than COST(s). 
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Many problems require that all the solutions satisfy 
some set of constraints which may be divided into two 

categories/ explicit and implicit. If a solution is 

represented by an n-tuple ( / X 2 / . . . / ) / then the explicit 

constraints are ones which determine what values the x^'s 

may take on. The implicit constraints restrict the ways in 
which the x^. 's in a solution relate to each other. All solu- 
tions which satisfy the explicit constraints are said to De- 
long to the possible solution space. The implicit con- 
straints then determine which solutions in the solution 
space satisfy the constraints of the problem. These solu- 
tions make up the 'feasible' solution space. 

For the graph colorability problem/ we will use the con- 
vention that a solution is expressed as an n-tuple 

( X. / X. / . . . / X ) of integers where l< = x.< = n/ and x. is the 
1 2 n 1 1 

color assigned to vertex i. The explicit constraints res- 
trict the x^. 's to integer values from I to n. Therefore/ 
the possible solution space has n'^ elements(n ways of pick- 
ing an integer from 1 to n for each of n positions). The 
implicit constraint in this problem is: 

(i / j) 6 E => X. J! X. 

1 J 

The size of the feasible solution space is therefore 
also a function of the number of edges in the graph and 
their arrangement. 
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The first step of branch and bound is to systematically 



divide the 


solution 


space into 


subset s . 


This 


may 


be 


represent ed 


by a tree 


organ i zat i on / 


where the 


root 


of 


the 



tree represents all feasible solutions. The children of a 
node N then represent the subsets into which N can be divid- 
ed. The leaf/ or terminal nodes are the particular solu- 
tions which make up the set of feasible solutions. 

For the graph colorability problem/ the following tree 
organization was used. The root node represents all feasi- 
ble n-tuples. Then/ an arbitrary vertex is chosen and the 
n-tuples are divided into n subsets such that subset(k) con- 
tains all the n-tuples for which x^ = k. This is then done 
repeatedly/ picking a new vertex at each level of the tree. 
For notation purposes/ a partial solution ( x^ / X 2 / . . . / Xj^ ) / 
k < n/ is an assignment of colors to the vertices 1/2/. ../k. 
This k-tuple is then used to represent all the elements of 
a subset. Using this notation/ a portion of the tree for 
the case n = ^ is shown in fig. 1. Note that the lower 
bound on the cost of a node is the number of different 
colors in the partial solution. 

A branching strategy is a rule for determining how a 
solution is to be divided into subsets. For this problem 
the branching function is the method of choosing the next 
vertex to color at each level of the tree. 
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( 1 ) 



( 2 ) 
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\ / \ 



/V 

/ \ 



( 1 , 2 , 2 , 1 ) ( 1 , 2 , 2 , 2 ) 




The tree structure for the Colorability Problem. 

Figure 1 



Once the organization is decided upon, the task of the 
problem solving method is to explore, or searcn this 'state 
space tree' until the optimal answer state is found. Since 
the time required by the search algorithm is a function of 
the'number of nodes exolored in the tree, the search stra- 
tegy is to explore the fewest nodes possible. One way to 
facilitate this is to 'prune' the tree as early in the 
search as oossible. Pruning takes place when a node is 
reached which violates the implicit constraints of the prob- 
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lem. Since no legal answer state could possibly be reached 
through that node^ none of its children need be explored. 
The size of the search space is thus reauced by the number 
of nodes below the pruned node. Obviouslyr the higher in 
the tree that pruning takes places the greater the reduction 
in the search space. 

The second way to reduce the search space is the essence 
of the branch and bound method. Each nooe has a cost asso- 
ciated with it» and since the cost function is non- 
decreasingf no searching need take place below any node 
whose cost is greater than the cost of some known solution. 
A 'search strategy' is a rule for choosing which unexplored 
node to explore next in the tree search. Several different 
basic search strategies exist for exploring a tree. 

A deoth first search proceeds down the tree from the 
root# exploring the leftmost children of each node first. 
At any point where pruning takes place# the search is backed 
up to the first node where it can again proceed down the 
tree in a different path. Top-to-bottom# 1 ef t-to-r i ght is a 
good description of this method. 

A breadth first search explores all the children of a 
node before oroceeding down the tree. Tnis might be called 
a 1 e f t - t o-r i gh t # top-to-bot tom search. 
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A *best-first' or 'least cost' search strategy exolores 
only the most promising path through the tree at any given 
time. A node is explored by assigning costs to all of its 
children and placing them on an unexplored list. The next 
node to be explored is the one on the unexplored list whose 
cost is currently the lowest. A priority gueue can be used 
to maintain an ordering of the unexplored nodes and at the 
point where the highest priority node is an answer state 
then the search can be terminated because that node 
represents the best or optimal solution to the problem. 
Nodes which violate the problem constraints may be assigned 
an arbitrarily high costr or merely not placed in the gueue» 
which effectively prunes the tree. This is the search stra- 
t egy employed in this research effort. 

This has been a general description of the branch and 
bound technigue. More specific details of the tree arrange- 
ments the search strategy# and the cost fu'nction used in 
this research will be described in Section III. 

C. PROBLEM complexity 

The complexity of a problem is said to be polynomial- 
time if an algorithm exists for which the number of elemen- 
tary or fundamental operations is limited by a polynomial of 
the length of an encoding of the problem. A problem is 
oolynomial-soace if an algorithm exists for which the amount 
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of elementary computer storage space needed for computation 
is never greater than some polynomial function of the encod- 
ing of the problem. An algorithm for a po 1 ynom i a 1 -t i me 
problem is said to be a 'good' algorithm and all problems 
for which a good algorithm exist belong together in a class 
called P. It may seem somewhat extravagant to group all 
polynomial-time algorithms into one class since polynomials 
can be quite large. However, no matter what the coeffi- 
cients are, in the limit as n gets large, every exponential 
function dominates every polynomial function. Further, 
though unexplained, experience has shown that for many of 
the problems encountered which have a po 1 ynom i a 1 - 1 i me algo- 
rithm, the solution is bounded by a polynomial of small de- 
gree. Some examples are ordered searching which is 0(ln n), 
sorting which is 0(n*ln(n)), and matrix multiplication, 
which is 0(n^"^^). 

Another group of problems are those whose best 'known' 
algorithm require greater than polynomial time. Some of 
these are; 1) Algorithms requiring subexponential time, such 

r 

as 0(e^^), 2) algorithms requiring exponential time, such as 
Ofe*^), and 3) algorithms requiring super-exponential time, 

en 

such as 0(e ). The disadvantage of non-po 1 ynom i a 1 -t i me al- 
gorithms is in the explosive growth of the maximum solution 
time as illustrated by figure 2, In this table the maximum 
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Comparison of several polynomial and exponential time 
comolexity funct i ons [Ref .51 
Figure 2 

Note that the only concern is with known algorithms and 
worst-case situations. No conclusions are made about ex- 
pected or average performance of an algorithm. 

Problems for which there exists an algorithm which can 
guess a solution and verify its correctness in polynomial- 
time make up the class of non-de t e rm i n i s t i c po 1 ynom i a 1 -t i me 
problemSf NP, For example# graph coloring is in NP because 
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an algorithm exists which can guess a coloring and check 
whether or not it is legal in po 1 y nom i a 1 -t i me , Such an al- 
gorithm needs 0(n) time to aroitrarily assign a color to 
each of the n vertices. Then it needs to consider each edge 
of the graoh and check whether the 2 vertices have the same 
color. This can be done in constant time for each edge and 
since there may be 0(n2) edges^ the entire verification 
takes 0 (n2 ) , 

Note in all of this discussion the word 'known'. There 
is not one single problem in the class NP for which it has 
been proved that the problem is not in PtRef. 6]. There are 
no po 1 y nom i a 1 -bounded algorithms known for many oroblems in 
NP/ but no lower bound which is greater than polynomial time 
has been proven for. any of these problems. Thus/ it is 
still an open question as to whether the class NP properly 
contains P or whether the two classes are equal. 

There is obviously much interest in trying to prove 
whether a given problem Q could or could not be solved in 
po 1 ynom i a 1 - 1 i me . For if one can prove that no polynomial- 
time algorithm for Q could Possibly exist then there is no 
point in expending the effort to find such an algorithm. 
Short of proving that a problem has no po 1 ynom i a 1 -t i me algo- 
rithm/ however/ there is some comfort in knowing that ones 
difficult problem is somehow related or at least as diffi- 
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cult as many other problems for which no one has found a 
po 1 ynom i a 1 *t i me algorithm. For this discussion some more 
definitions are needed. A 'problem' is some general ques- 
tion to be answered which usually has one or more unspeci- 
fied parameters and some specified properties which the 
answer is required to satisfy. An 'instance' of a problem 
is obtained when the parameters are specified. For examolef 
an instance of the colorability problem is obtained by 
specifying the vertices and edges of the graph and the 
number of colors for which a coloring is desired. 

A problem Q is said to be 'reducible' or 'transform- 
able' to R if an instance of Q can be transformed in poly- 
nomial time to an instance of R. In other wordSr if the 
answer to problem Q is 'yes' if and only if the answer to R 
is 'yes' and if the input parameters to problem R can be 
determined in polynomial-time given the input parameters to 
Qf then Q is reaucible to R. The following schematic illus- 
trates this definition; 



Input 
to Q 



Po 1 ynomi a 1 


Input 


A 1 gor i t hm 


t i me 


to R 




t ransform 




for R 


to R 







Output 
for R 
and Q 



> 



Transformation of Problem Q to Problem R 

Figure 3 
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Clearly^ if Q is reducible to R and R is in the class P# 
then Q is also in P. Note that nothing is said about the 
existence of an algorithm for R» only that 'if* an algorithm 
for R exists which runs in oolynomial-time/ then one also 
exists for Q. 

If every problem in the class NP is reducible to some 
problem Q then Q is said to be NP-hard. In other words^ Q 
is at least as hard as every problem in NP, Furtherf if Q 
is also in the class NP then Q is NP-complete. The signifi- 
cance of a problem being NP-complete is that if a 
po 1 ynom i a 1 - 1 i me algorithm is ever found for any NP-complete 
problem then it will be known that a oo 1 ynom i a 1 -t i me algo- 
rithm exists for all NP problems and that indeed NP = P. 
Further/ if any problem in NP is ever shown to require 
super-oo 1 ynom i a 1 time complexity then every NP-complete 
proolem must also require greater than polynomial time. 

To show that a particular problem Q is NP-complete/ 
then/ one must show that it is indeed in NP and then either 
show that every NP problem could be reduced to Q or that 
some already proven NP-complete problem is reducible to Q, 
In 1971/ Steven CooktRef, 71 laid the ground work for the 
current theory of NP-comp 1 et eness in a paper entitled "The 
Complexity of theorem proving procedures," One of his most 
significant results was the proof that the 'Satisfiability' 
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problem from boolean logic is NP-completer and it holds the 
distinction of being the first NP-complete problem. Refer- 
ence 5 contains an excellent description of the Satisfiabil- 
ity problem as well as a proof of Cook's theorem. Cook also 
proved that a variation of the Satisfiability problem^ 
called 3-Satisfiability/ was NP-complete by transforming Sa- 
tisfiability into 3-Satisfiability. Then/ in 1972/ 
KarplRef. 8] proved that many decision problems were also 
NP-complete. One of these proofs was that the 3- 
Satisfiability problem was reducible to the colorability 
problem. Thus the colorability problem is also NP-comolete. 

Inasmuch as most theorists are in agreement that the 
NP-complete problems are probably super-polynomial/ much 
work has been done in recent years to find relaxations to 
the problems which allow them to be solved in polynomial 
time. As the graph colorability problem has many apolica- 
tions in real-life problems/ many methods of finding approx- 
imate solutions have been tried. In many cases where a 
graph can be used to model a real-life problem/ the absolute 
minimum number of colors may not be needed/ but rather/ some 
close approximation might be satisfactory. Unfortunately/ 
even this has proved to be a difficult problem. Many algo- 
rithms have been dev i sed [Re f s , 2/3/^U to find approximate 
colorings/ but for every one/ it is possible to construct a 
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graph such that the approximation algorithm will find a 
coloring which uses many more colors than the chromatic 
number , Though the algorithm may indeed run in polynomial 
time» the closeness of the approximation to the actual 
chromatic number of the graph may be suspect. In fact» in 
1976# Garey and JohnsonCRef. 1^1 proved the following; 

Given graoh G with chromatic number X(G)# and 
po 1 ynom i a 1 "t i me approximation algorithm A» which computes a 
coloring of G using A(G) colors# there can be no A such that 
for all G : 



A(G) 

< r for r < 2, 

X(G) 

In other words# any approximation algorithm will always 
have some input graph for which it can find no coloring 
which uses less than twice as many colors as actually need- 
ed. 

The primary reference for the material in this section 
is Garey and JohnsonlRef, 51. It is highly recommended for 
the reader who is interested in other NP-complete problems 
and further discussion on this topic. 
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THE ALGORITHMS 



A. GRAPH GENERATION 

One of the earliest problentis faced in this research was 
that of deciding what characteristics were desirable in the 
graphs to be colored and how to generate them. As the goal 
of the research was to contrast the efficiencies of various 
heuristics for optimal colorings# the graphs used needed to 
be# in some sense# typical# or representative of graphs 
which occur in real-life problems. Though the ouestion of 
what is a tyoical graph was left largely unanswered# three 
possibilities for graph generation were considered. 

One method considered was to use hand generated graphs 
with certain built-in specific features. For example# 
graphs could be designed with particular distributions of 
degree of the vertices. Another method would be to generate 
graphs with a given chromatic number but with varying 
numbers of vertices or edges. The major drawbacks to this 
method were the necessarily small sample of graphs and the 
failure of the graphs to be typical in any sense. This 
method was therefore not used for any of the experiments 
done but was used extensively for program testing and debug- 
ging. 

Another possibility# guickly discarded# was to survey 
some of the disciplines where graph colorability models the 
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Dpoblem and to use some actual graphs which occur in real- 
life. Though this might generate some interesting results 
for those specific applications^ it would also lead to 
results which were heavily flavored by the problem being 
considered. rthat was needed was a selection of graphs which 
were# in some senses representat i ve of the set of all 
graphs . 

The most promising method of graph generation was to use 
the computer and to create a large population of random 
graphs on which to test the various heuristics. There is 
some benefit to using random graphs since one generally 
feels comfortable that these random graphs are fairly 
representative of the set of all graphs. 

The primary drawback to random generation of graphs is 
that the generator will tend to generate more of the types 
of graph which have many isomorphisms than the type with few 
isomorphisms. A graph G is isomorpnic to graph H if the 
vertex numbers of 6 can be permuted in such a way as to make 
the resulting adjacency matrix of G identical to the adja- 
cency matrix of H, Fig H shows 2 of the isomorphic graphs 
for each of the only 2 distinct graphs possible for ‘i ver- 
tices and edges. 

Since graph coloring algorithms can be influenced by the 
different types of graphs^ a random generator may bias the 
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investigation of the relative performances of different 
coloring programs. The generation technique could be amend- 
ed by somehow reducing the probability of using a particular 
graph by dividing the probability by the number of the 
graph’s isomorphisms. But determining the number of iso- 
morphic configurations of a given graph is not a straight- 
forward problem and beyond the scope of this research. 





Examples of the isomorphs of 2 distinct graphs. 

Figure 4 

Another technique might be to record each graph and to 
test each new one generated against all previous to deter- 
mine whether it is an isomorphism of an already used graph. 
But the problem of just testing two graphs against each oth- 
er is probably NP-compl ete( though considered an open ques- 
tion at this time) CRef. 51. 

In order to get on with the research goals at hand# it 
was decided to accept a random graph generator which gen- 
erates graph type G more than type H if G has more isomor- 
phisms than H. This is an admitted shortfall in this 
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research effort and it is hooed that the algorithms experi- 
mented with are not discriminated against by the shortcoming 
of this random graph generator* 

Inasmuch as unconnected graphs are no more difficult to 
color than their largest connected componentf it was decided 
to run the initial tests on only connected graphs. Two 
methods of generating connected graphs were considered: 1) 
To randomly add edges until the graph was connected and 2) 
to create a spanning tree of n vertices and n-1 edges(the 
minimum to connect the n vertices) and then add additional 
edges at random* Since method one would require testing for 
connectivity after each edge and also would fail to generate 
very many sparse graohsCfound often in real oroblems) method 
2 was chosen and implimented using the following algorithm: 



Algorithm GENGRAPH 

1. Let VC be a set of connected vertices^ initially empty 
Let VU be a set of unconnected vertices# initially 

VU = {l#2#..*#n> 

2. Choose a random vertex a 

3. Move a from VU to VC 

• While VU is not empty 

choose a random vertex a from VC 

choose a random vertex b from VU 

record an edge from a to b 
move b from VU to VC 

5. While Number edges < Number desired 
choose random vertex a 

choose random vertex b 

While edge (a»b) exists or i f a = b 

increment a until a = n then increment b and set a = 1 
if a = n and b = n then set a - 1 and b = 1 
record edge (a#b) 

End 
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At a later point in the research# it was decided that 



for truly valid test results# the graoh generator must gen- 
erate purely random graphs with no requirement for connec- 
tivity. It was discovered that the connectivity requirement 
changed the characteri st ics of the sparse graphs and non- 
represent at i ve graphs were being used. Therefore# a much 
simplified generator was used which ignored steps I through 
^ of 6ENGRAPH and merely added random edges until the re- 
quired number was oresent. Another algorithm called DIVIDE 
(to be described in a later subsection) was then called to 
divide the graph into its connected components. 

The random number generator used is identical to Grogono 
(Ref. 91 with a modification to return integers in any range 
desired. The listing for both GENGRAPH and the Random 
number generator is enclosed in appendix A. 
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B. GRAPH DIVIDING 



It has already been shown that an unconnected graph may 
be colored by dividing it into its connected components and 
separately coloring each component. A further simplifies/ 
tion can be made# however/ if even the connected components 
are divided at their articulation points/ if any exist. Re- 
call that an articulation point is a vertex which is in 
every path from some vertex x to some other vertex y. If the 
articulation point is divided into 2 vertices such that the 
oath from x to y is cut/ then the graph is divided into 
disconnected subgraphs. For illustration purposes consider 
figs. 5 and 6. 




A Graph with articulation points. 
Figure 5 
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The articulation points are vertices d and f 



The 



biconnected components are shown in fig, 6. 




A Graph Divided at its articulation points. 

Figure o 

In communications and transportation problems the iden- 
tification of articulation points is very important because 
these are the 'choke* ooints» or places where the network 
becomes most vulnerable to failures. Since dividing a graph 
at its articulation points results in biconnected components 
which can be colored separately^ it is also important to 
discover the articulation points in the colorability prob- 
lem, The number of colors needed will be exactly the same 
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as the number needed to color the component which uses the 



most colors. If the actual color assignments are needed/ 
then the graoh can be 'pieced’ back together one component 
at a time. As each component is added/ the color classes 
are permuted so that the articulation vertex is in the same 
class in both components of the graoh. 

An algorithm to find the articulation points and the 
biconnected components of a graph utilizes a depth first 
search of the vertices in the graph and mai<es use of the 
fact that V is an articulation point if and only if there 
are two vertices x and y such that every oath from x to y 
includes v. 

The idea of the depth first search is to visit all the 
vertices of the graoh by starting at some arbitrary start 
vertex s and recursively visiting every other vertex. Given 
that the search is at some vertex a, it follows an edge 
(a/b) to vertex b. If b has already been visited/ the 
search returns to a and tries another edge. If b has not 
been visited/ then the same method is applied recursively to 
b. At the point where all the edges incident on a have been 
examined/ the search backs up along the edge that took it to 
a. The search terminates at the point where it tries to 
back uo from the start vertex. If the graph was not con- 
nected/ then some new start vertex is chosen from the un- 
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visited ones and the procedure repeated until all vertices 
have been visited. The edges which take the search from one 
vertex to a vertex which has previously been visited are 
called 'back edges.' Edges which go to an unvisited vertex 
are called 'forward edges.' These definitions will be need- 
ed for the algorithm description. 

Figure 7 shows a graph with 5 vertices and 7 edges. If 
the depth first search is started at vertex a then the solid 
arrows represent a possible set of forward edges and the 
dashed arrows are the back edges. The labels on each edge 
indicate the order in which they are traversed. 




A Depth first search of a Graph. 

Figure 7 

The basic idea of using a depth first search to find the 
articulation points of a graph can be seen by studying fig. 
8. This schematic represents a graoh with 3 piconnected 
components labeled and C which are joined at articula- 

tion vertices v and w. S is an arbitrary start vertex in A. 
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A Schematic of a graoh with 3 biconnected components and 

2 articulation points. 

Figure 8 

If the depth first search is started at vertex S in Af 
it will eventually enter into C through vertex v. By the 
depth first nature of the search/ all the edges in C will be 
traversed before the search backs up through v. If the 
edges are placed on a stack as they are traversed/ then wren 
the search gets back to v/ all the edges in C will be on top 
of the stack. 

In order to recognize the articulation points/ note that 
no vertex in C will have a back edge to any vertex which was 
visited prior to v, So/ if the vertices are numbered as 
they are visited and each vertex is tagged with an index 
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equal to the smallest number of any vertex which can be 
reached through any number of forward edges and one back 
edge» then the articulation points can be found. When the 
search backs up from b to a along edge (afb)/ if the index 
tagged to b is not less than the number assigned ar then a 
is an articulation point(or the root of the tree). In other 
words/ no vertex below a in the tree has a back edge to a 
vertex visited before a. The following algorithm by 
BaaselRef. 6) was used to find both the connected components 
and the biconnected components of a graph: 

Algorithm DIVIDE 



SV = stack of vertices 

SE = stack of edges 

top = top element of SV 

Number = array to number the vertices in the order visited 
Back = array to record the last vertex reachable 
through a back edge. 

Algorithm: 

1. initialize Number(i) = 0? 

2. choose arbitrary vertex S 

3. Number(S) = 1 

y. Num =2 // next number 

5. stack S on SV (top = S) 

6. while there exists an unprocessed edge from top to v do 

stack edge (too/v) on SE 

if number(v) > 0 then back(too) = m i n (bac k ( top) / bac k ( v ) ) 
else do 

number(b) = num 
num + + 

back(v) = number(v) 

stack V on Sv 

end 

end 
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7. if there is more than one vertex on SV 



let V = second from top 
if back(top) >= Number(v) then 

let edge (a^b) be the top edge on S£ 

while Number(a) >= Number(v) and Number(b) >= 
Number ( v ) 

POP S£ 

else Back(v) = m i n (Bac k ( v ) » Bac k ( t op) ) 

POP SV 
goto 6 £N0 

The code used to impliment this algorithm is enclosed in 
appendi x A , 

C. GRAPH COLORING ALGORITHM 

The essential idea of the basic coloring algorithm is to 
build the search tree until the optimal coloring is found. 
£ach node in the tree represents the placement of a vertex 
into a color class and the oath from the root to the node 
represents a partial coloring of the graph. The cost of a 
node is the number of different colors used for that partial 
coloring. As the nodes of the tree are generated^ bv itera- 
tively expanding from the start nodcf they are placed in a 
priority gueue such that the first node in the gueue is the 
node which represents the least numoer of colors used and is 
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deeper in the tree than any other node with the same number 
o^ colors. A node is expanded by generating all of its pos- 
sible children and the node to be expanded next is always 
the first node in the queue. The stopping point for the al- 
gorithm is thus the point at which the highest priority node 
is a full coloring of the graph. At that point no other 
node in the tree could possibly be expanded in any way which 
results in a coloring using fewer colors. The chromatic 
number of the graph is then known to be the number of colors 
used in that coloring. Note that if all possible colorings 
are desired which use the optimal number of colors/ then the 
algorithm need only save the chromatic number and to contin- 
ue expanding nodes until the cost of the highest priority 
node is greater than the chromatic number. At each point 
where a node is generated which is a full coloring/ it can 
be printed out and all optimal coloring assignments will 
thus be generated. 

There are 2 primary data structures used in this algo- 
rithm. These are the TREE and QUEUE data structures. The 
elements of TREE are trees made up of nodes and the links 
between them which organize the nodes into a tree structure. 
A node N represents a color assignment to one vertex of the 
graph. The path from N to the root of the tree represents a 
partial coloring of the graph and the cost assigned to node 
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N is the number of colors used in that partial coloring. 
Each node has a pointer to its parent in the tree in order 
to maintain the links between nodes. There are 2 operations 
defined for TREE, The operation CREATE generates a tree 
with only one node called the root. The operation EXPAND 
takes as input a tree T and a particular leaf node N of T 
and returns a tree T which is the original tree T plus all 
the possible children of node N, 

The nodes of the tree are imolimented using an array of 
records where each record contains the following 5 fields: 

1) the depths or level in the tree at which it is located^ 

2) the cost assigned to the node/ 3) the vertex number which 
was colored to create the node/ 0) the color assigned to 
that vertex/ and 5) a pointer to the node’s parent in the 
tree. The operations on the TREE are implimented in the al- 
gorithm COLORALL and EXPAND as follows: 



COLORALL (Color all vertices of the graph) 

1 . Get a node N 

2. Get number of first vertex - v . 

3. Let N be the root node 
N. vertex = v 

N.cost = 1 
N, depth = 1 
N . CO 1 or = 1 
N, parent oointer = NULL 
y. While N, depth < n do 

expand N (see algorithm EXPAND below) 

N = first node in priority queue 
5. Print out s t a t i s t i c s / co 1 or i ng / and chromatic number. 



The algorithm to expand a node is as follows: 



EXPAND(N) (N is the node to be expanded) 

1. Choose a vertex v to be added to the partial coloring, 
(described in subsection D.) 

2. let C = { 1 , 2» . . . / N .cost > be the set of colors used 

in the path to node N. 

3. Trace the parent pointers from N to the root. For 
each node M in the oath to the root; if an edge exists 
between vertex v and M. vertex then remove color M, color 

f rom C . 

For each color c- left in C create a new node N. and 



END 



put it in the priority queue. 
N. .vertex = v 



,1 



= N.cost 



N. .cost 

n]. depth = N. depth + 
nI .col or = c. 



.1 



r = 



.parpt I 

Aad one new color and create one more node N 

N..vertex=v 

N'^.cost = N.cost +1 

N'^. depth = N. depth +1 

N^ .CO 1 or = N.cost t I 

N'? .parot r = N 
0 



In terms of graph coloring, step 4 creates a child 
node for each of the ways in which the new vertex can be ad- 
ded to one of the existing color classes. Step 5 generates 
a node which represents putting the new vertex in a new 
color class by itself. 

The QUEUE data structure is used to organize the leaf 
nodes of a tree according to their priority. There are 2 
components which determine the priority of a node. The 
first is, of course, its cost as defined earlier. The 
second component is used to determine the highest priority 



between 2 nodes with the same cost 



When more than one node 



is available which have the same cost then the node which is 
deepest in the tree has a higher priority since it is closer 
to a complete coloring. There are two operations defined 
for the QUEUE cal led ADDNOOE and REMOVETOP, AOONODE takes a 
queue Q and a node N and returns a queue Q with N in its 
proper priority location, REMOVETOP takes a queue Q and re- 
turns the node N with highest priority. 

The priority queue is implimented as a heap# as 
described in ref, 10, The nodes are organized into a full 
binary tree such that the root is the node with highest 
priority and every node in the tree has a higher priority 
than either of its children. When adding a node to the 
tree» it is put in the first available location and then 
•sifted up* until its parent has a higher priority. When 
the root node is removed# the last node is moved to the root 
position# then ’sifted down* until both its children have a 
lower priority. The binary tree is maintained by an array 
of pointers q(l where q(l) Points to the node at the root of 
the binary tree and the children of q(i) are located at 
q(2i) and q(2i + 1), 

In any tree search algorithm it is important to ensure 
that isomorphic paths are not generated. For the problem of 
graph colorability# two colorings are isomorphic if the ver- 



tices of the graph are partitioned into color classes in the 
same way in both colorings. Therefore* if one coloring uses 
k color classes* then there are kl possible isomorohs of 
that coloring. In order to prevent the tree from containing 
all of these unneeded paths* one must either check for pos- 
sible isomorphs with the addition of each new node or 
prevent isomorphic oath generation from occuring by the 
design of the algorithm. 

There are two essential ingredients of algorithms 
COLORALL and EXPAND which ensure no isomorphs: 1) When a 
node N is expanded* at most I new color is added to the 
available color set* and 2) In expanding a node N* only I 
vertex of the uncolored vertices can be added to the partial 
coloring and that vertex must be the same for all the chil- 
dren of N. 

THEOREM; Any tree search graph coloring algorithm with the 
above 2 properties will not generate isomorphic colorings. 

Proof: Assume the algorithm generates a tree with 2 iso- 
morphic paths P and Q. Then* by the definition of isomorph- 
ic colorings* for each color class in P there exists a color 
class in Q with exactly the same elements. Let vertex v be 
the first vertex colored in P and Q such that the color of v 
in P is different from the color of v in Q. Call the class 
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p 

containing v in path P/ C!j^ and the class containing v in Qf 
q 

C ,• . If the subscript indicates the color/ then i X j. 

Since only one new color could be added when v was colored/ 

P Q 

then at least one of or Cj contains a vertex w which was 
colored prior to v, without loss of generality/ let that 
class be C^. In other words/ v and w are in the same color 
class in path P. But/ since v was colored differently in 
oath Q and w is the same in both P and Q/ there can be no 
color class in Q containing both v and w. Therefore/ P and 
Q are not isomorphic colorings. Figure 9 illustrates this 
proof . 
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Two Paths P and Q which are not isomorphic 

Figure 9 



D. THE VERTEX ORDERING FUNCTIONS 



The purpose of this research was to investigate dif- 
ferent branching functions in a branch and bound approach to 
solving the chromatic number problem. For the particular 
branch and bound algorithm used# the branching function is 
the order in which the vertices are selected to be added to 
the partial coloring at each new level of the tree. In this 
research# five different methods of ordering the vertices 
were investigated. Four of the methods were static order- 
ings in which the vertices were preoraered and then the 
graph coloring algorithm was called to search the entire 
tree. In every path through the tree in these methods the 
order of choosing the next vertex was constant. In the 
fifth method# the vertex chosen to be colored next was 
dependent on the particular color assignments made prior to 
the node to be expanded. This method of dynamic ordering 
results in a larger time expenditure at each node since the 
order must be recomputed at each expansion of the tree. 

Note that the purpose of any of these orderings is to 
make the search tree as small as possible# thus requiring 
fewer nodes and fewer expansions. If the amount of work 
done to expand a node (create its children) is constant for 
each method# then the resulting size of the search tree is a 
direct measure of the benefit of any procedure. Since the 
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height of the tree will always be the number of vertices in 
the graoh, the only way todecrease the number of nodes is 
to decrease the branching factor (number of children of each 
node). The optimal branching factor is# of course# 1. With 
the algorithm COLORALL# this is reached when each vertex to 
be colored requires a new color. It is approached most 
closely when the vertices are ordered such that the next 
vertex is always the one with the fewest colors available to 
be colored. Three of the static orderings attempt to accom- 
plish this ordering through the use of some heuristic func- 
tion to order the vertices. The dynanvic ordering actually 
determines# at each node# which of the uncolored vertices 
has the fewest colors available. 

Obviously# any graph which is complete(a clique)# will 
be colored with a branching factor of one# no matter what 
order is used. The same applies to a graph which is just a 
forest of vertices(no edges). The difficulty lies in the 
midrange. Another way to look at it is that coloring is 
easiest when there are few res t r i c t i ons (edges ) or very many 
restrictions to the colors which can be used. Each of the 
ordering methods# therefore# really attempts to find the 
vertices with the most restrictions and color them first. 
Thus# the tree has the smallest branching factor at the top 
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and more of the tree is pruned off resulting in fewer nodes 
to search. 

Following is a description of each ordering. The list- 
ing for all orderings is in the program ORDERNODE in appen- 
dix A, The case numbers in the switch statement correspond 
to the numbers shown here, 

1 , Random 

This method of selecting vertices was used as a 
baseline against which to compare any benefit of the other 
methods. Inasmuch as the graph was generated randomly^ no 
real work needed to be done. The vertices were merely 
chosen in the order of increasing vertex number, 

2 , Degree 

This method# used in the Vfe 1 sh-Powe 1 1 approx i mat i on 
a 1 gor i t hm IRef , 111# orders the vertices by decreasing values 
of the degree of the vertices. Recall that the degree of a 
vertex is the number of edges incident upon it and is easily 
calculated by summing the elements of the rows of the adja- 
cency matrix. The idea behind this method is that the ver- 
tex with the highest degree has more restrictions on it when 
choosing a color so it should be colored first. That color 
used is then removed from the colors available for every 
vertex connected to it. 



^<7 



3 . Dvec 3 



OVEC 3 is the name given by this author to a method 
devised by M.R.Wi 1 1 i ams IRef . 4]. In the prior method of 
ordering the vertices by degree » the underlying premise is 
that if a vertex i is connected to p other vertices^ then it 
will be more difficult to color than vertex j which is only 
connected to q vertices (p < q) , What that method fails to 
consider is the degreer or difficulty^ of the vertices to 
which i and j are connected. For example^ if vertex j is in 
a clique of size 4 and i is the center of a star with 5 
points/ as shown in fig. 10/ then certainly j should be 
colored before i/ as well as all the other vertices in the 
c 1 i que . 




Example of Shortcoming of Degree method 

Figure 10 



Williams recognized that the vector of degrees of the ver- 
tices was the first step towards the development of the dom- 
inant eigenvector of the adjacency matrix. Recall that 

d. = Z a. . 

1 ; 1 J 
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In the standard technique/ d^- is the eigenvector 
corresponding to the largest eigenvalue of matrix A and is 
found by iteration using: 



n + 1 



= I 

J 



'<1 



If d^ is the vector of degrees/ then this scheme 
does take into account the degrees of the vertices to which 
a vertex is connected. And/ in fact/ the number of itera- 
tions used determines how far away from vertex i to look in 
calculating the relative difficulty of coloring vertex i, 

Williams also determined that the relative order of 
after a few iterations. This fact was borne out by results 
of this research. Although improvement was seen out to 2 or 
3 iterations/ little or no changes in the size of the search 
tree was found after 3 iterations. The only effect was the 
increased computation to calculate the vectors. For this 
reason/ 3 iterations was chosen as a good number to use to 
compare this method against the others. 

4 . Ad i acencv 

The rationale behind this ordering method is that 
the vertex which is adjacent to the most previously colored 
vertices will probably have the fewest colors available. 



This is not necessarily true» however^ as a vertex may be 
adjacent to many vertices which are all in the same color 
class» thus reducing its available color space by only one 
color. 

The algorithm to determine which vertex is next# 
considers the rows in the adjacency matrix for all the pre- 
viously colored vertices# and adds the columns. The highest 
column total then indicates which vertex is to be colored 
next. Previously colored vertices can be eliminated by sub- 
tracting an arbitrarily high number from the running column 
total# for each vertex# when it is chosen. In the case of 
ties# the vertex of highest degree is chosen. 

5 . Least colors available (LCA) 

In all the previous static orderings# an assumption 
was made that the best vertex to color next# at any point in 
the tree# was independent of the specific color asignments 
made up to that point. We were simply trying to find the 
vertex that was 'most probably* a better choice than the 
others. In this dynamic ordering# when expanding each node# 
we actually choose the vertex which will result in the 
fewest children being created# thus minimizing the branching 
factor. Since a child node is created for each available 
color which the next vertex can be colored# we must consider 
all uncolored vertices in the path to the node being expand- 
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ed and determine which vertex has the fewest colors avail- 
able. To do this note that ideally we would like to choose 
a vertex which is adjacent to at least one vertex in each 
color class used so far. Thus a new class would have to be 
created and there would be only one child node created. 

Short of this ideal situationr we would like to choose the 
vertex which is adjacent to at least one vertex in most of 
the color classes so far created. 

This is accomolished by creating a new matrix CC 
similar to the adjacency matrix such that there is a row in 
CC for each color class generated and a column for each ver- 
tex. Then we define 

CC . . = j 1 if (Vfj) 6 E for some vertex v in color class i. 

^ I 0 otherwi se . 

i.e.f CC . . = 1 iff vertex j cannot be colored color i. 



The matrix CC is easily created from the adjacency 

matrix A as follows: Let cc . be row i in CC/ and a be row i 

1 i 

in A, both represented as bit strings. 

initialize all cc . . = 0 

for i from 1 to Nil3|-C0L0R-CL ASSES 

for each vertex j in color class i 
c c = c c V a 

The next step is to generate a vector called SUM such that 



SUM .= 
J 




cc 



1 J 



if j 



i f j is 



i s uncol ored 
col ored 
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Then the value of j for which SUM is a maximum is 
the vertex which has the fewest color classes available to 
which it can be added. Again/ for consistency/ ties are 
broKen by choosing the vertex of highest degree. 

E. CONTROLLING ALGORITHM 

Inasmuch as the ourpose of this research was to experi- 
ment with various orderings of the vertices/ the flow of 
control was fairly st ra i ght f orward. After only a few trial 
cases it was discovered that the primary limitation to the 
size of the graph which could be colored was computer 
storage space/ not time. As the graphs got larger/ more 
nodes were required in the search tree until the available 
memory was used up at about ^000 nodes. Although unsigned 
variables were used in the node structures and space was 
conserved wherever possible/ this limitation could not be 
significantly changed. For even the best orderings/ this 
meant graphs of about 30 vertices or less. For empirical 
data then/ graphs were generated with from 10 to 30 vertices 
and edges from BOX to 80% of the maximum possible. For each 
graph size/ 100 graphs were generated (some of which were 
probably identical and many which were at least isomorphic) 
and each ordering method was used to color the graph. The 
graph was then divided at its articulation points and each 
bi connected subgraph was colored using each ordering func* 
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tion. Data about the size and shape of the search tree was 
collected and tabulated for all the above samples 
( append i x 0 ) . 

Note that 100 trials of each graph size is a relatively 
small sample^ but was a necessary restriction to keep the 
entire program execution time to within a few hours. 

The following algorithm shows the setup used for obtain- 
ing the experimental data. 



algorithm MAIN 



initial i ze 

For n from minsize to maxsize (n is number of vertices) 
begi n 

compute maxedges = (n)(n-l)/2 

For fullness from 30% to 70% of maxedges 

begi n 

compute number of edges e = fullness * maxedges 
For trials from 1 to 100 
begi n 

generate random graph (n vertices^ e edges) 

For each ordering method 
beg i n 

color the graph 

record data (expansions# nodes# branching factor) 
divide the graoh 

For each biconnected component of the graph 
begin 

color the component 
record component data 

end 

sum the component data 

end 

end 

compute averages of all data collected 
print out data 

end 

END 
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The function MAIN in appendix A is the controlling algo- 
rithm for the entire process, 

F, GENERAL COMMENTS ON ALGORITHM DEVELOPMENT AND THE 

PROGRAMMING LANGUAGE ’C’ 

In all the programming done during this research# ef- 
forts were made to save time and storage space whenever 
feasible. The first attempt to design a coloring algorithm 
used a backtracking method which is common for tree 
searches. However# the sire of the search tree was so large 
that the storage required to stack the variables used in 
each recursive call to the backtracking function was exorbi- 
tant, For that reason# an iterative scheme was chosen and 
i mol i mented. 

Further savings were made by using bit strings and 1 bit 
variables whenever possible. Since the maximum graph sire 
was to be less than 32 vertices# the long i nt eger variables 
of 32 bits in 'C' were ideally suited for the adjacency ma- 
trix and many other needed variables. Also 'C* has a good 
assortment of bit manipulation operators including multiple 
shifts# or# and# and exclusive or. Using these led to time 
savings in many cases. 

Space was conserved by using the unsigned variable type 
in arrays and records. With this convention# the length of 
the storage cells reserved were only as long as needed. For 
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examp 1 e t 



the vertex numbers from 0 to 31 could be stored in 



only 5 bits. Another practice^ not necessarily considered 
good programming technique was the fairly extensive use of 
global variables. Though risky^ it did save the time re- 
quired to pass parameters and the space to store local vari- 
ables. 

Overall/ this author feels that 'C' is a very appropri- 
ate language for research of this type. 
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IV. RESULTS 



A. THE PARAMETERS 

In attempting to determine the relative 'goodness* of 
several computer programs^ one must first find some measur* 
aPle parameters to be minimized or maximized to indicate the 
best program. For approximation algorithms^ one is very 
often not interested in how long a program takes as long as 
it is some linear function of the problem size or possibly 
even O(nlgn). What one measures is the closeness of the ap- 
proximation to the actual optimal solution. 

For this research/ the solution to the graph coloring 
problem was to be the exact chromatic number in every case. 
The coloring algorithm was to find the chromatic number of a 
given graph and the 'goodness' of the vertex ordering 
methods was to be determined by how fast the algorithm ran 
and how much space was used. The specific questions to be 
answered were! 1) Is there a best-time ordering for all" 
graphs?/ 3) Is there a best-space ordering for all graphs?/ 
3) Are they the same?/ What tyoe/ or size/ of graphs 
favor which ordering if there is a difference?/ and 5) Is 
one ordering better on the average/ but another better for 
worst-case situations. 

The results of the many experiments which relate to 
these questions is tabulated in appendix B. The discussion 
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here is of a general nature and will summarize the results. 
There were basically two separate categories of measurements 
made - those related to time and those related to space. 
Though they are closely related in this problem/ the parame- 
ters measured are distinct and will be discussed separately. 
1 . T i me 

In determining the difference in time-complexity of 
the 5 coloring methods/ there are 2 modules of the program 
to consider - OROERNOOE and C0L0RALL(of which EXPAND is a 
part). These functions are the only two whose time- 
complexity is a function of the vertex ordering as well as 
the graph size. 

For the ^ static orderings the ti me -comp 1 ex i t y is 

easy to determine. For consistency/ each ordering merely 

filled an array called next/ where next(i) was the vertex to 

follow i in the ordering. Random is therefore clearly 0(n) 

since the vertices are ordered by vertex number. One 

traversal through the array next is all that is needed to 

2 

fill it. DEGREE is 0(n ) since the adjacency matrix has to 
be tr.aversed to find the degree of each vertex/ then a sim- 
ole sorting of the vertices is O(nlgn) so the overall com- 

? O 

plexity is 0(n ), DVEC 3 is 0(n ) if the number of itera- 
tions is kept cons t an t ( sue h as 5 in this case). It would be 
conceivable to iterate n times for the best possible order- 
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3 

i ng in which case DVEC n would be 0(n ), ADJACENCY is also 

2 

0(n ) because a running total was kept of the relative oosi- 
tion of the unselected vertices as each new one was deter- 
mined. Thus the adjacency matrix was addressed in only one 
row for each new vertex. For the dynamic ordering LCA^ no 
work was done in the function ORDERNOOE but rather in the 
function PICKNODE which was called every time a node was to 
be expanded. PICKNODE itself has a time complexity of n^ 
since it takes n^ operations to create the n x n matrix of 
color classes CC and each element in CC requires 0(n) to 
compute. But the bigger concern is that the time-complexity 
of the branching function EXPAND is n^ times as great for 
LCA than for the static orderings. 

For the static orderings the time-complexity of gen- 
erating the search tree has 2 components - the time to ex- 
pand a node and the number of expansions needed to find a 
solution. For LCA# a third component# the time to find the 
next vertex# must also be considered. 

Once the next vertex is found# the time required by 
the function EXPAND is 0(n), The only iterative part of 
that function is following the path to the root and then 
generating a child node for each possible coloring. Since 
the depth of the tree can be up to n-1 and the number of 
children can be n# this function is 0(n), 
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So farf all the time complexities have been polyno- 
mial. If this were also true of the number of expansions# 
then one could claim that the overall algorithm would run in 
pol ynomi a 1 -t i me since the product of any number of polynomi- 
als is also a polynomial. Unfortunately# that is not the 
case for worst-case situations. Consider for a moment the 
size of the full search tree generated by this algorithm if 
no pruning is allowed. Refer to figure 11 which shows the 
first four levels of the search tree. The notation x#y in 
each node of the tree indicates that the node represents an 
assignment of color y to vertex x. Note that at the root 
there is 1 node followed by 2 nodes at level 2. Then at 
level three there are 5 nodes# 2 children of one in level 2 
and 3 for the other. To determine the number of nodes at 
level k# recall that one color class is added for each level 
in the tree. Also# a node at level k represents one of the 
ways of partitioning the vertices into k or fewer color 
classes. Therefore# the number of nodes at level k is the 
number of ways to partition a set of n elements into k or 
fewer classes. 

The Stirling numbers S(n#m) represent the number 
of ways to partition n elements into m distinct 
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THE FIRST FOUR LEVELS OF THE UNPRUNED SEARCH TREE 




Figure 11 
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classeslRef. 12] and are generated by the following recur- 
rance relation: 



i/n) = 1 1 

1,1) = 1 ) 



for all n 



S(n+l,k) = S(n,k-1) + k * S(n,k) 
where S(n 
and S(n 

To understand this relation, consider the following: 
Given that n elements can be partitioned into k classes in 
S(n,k) ways, there are two ways to add the (n+l)st element. 
First, it can be added to any one of the k classes for each 
way of partitioning n elements. The number of ways of doing 
this is given by k * S(n,k). Secondly, the n elements can 
be partitioned into k-1 classes and then the (n+l)st element 
can be a single element in a new class. This is the reason 
for the S(n,k-1) term. 

But the number of nodes at level k is the sum of the 
Stirling numbers S(n,m) for m from 1 to k. These numbers 
are called the Bell numbers B after E , T . Be 1 1 IRe f . 13). B 

n 

is the number of partitions of a set with n elements into 
any number of distinct classes up to n. That is 

= Z S(n,k) 
k 

To see how fast these numbers grow, the first 8 are: 

1 , 2 , 5 , 15 , 52 , 203 , 877 ,aiao 



In other words, at level 8 in the search tree there 
would be 41^10 nodes and the total number of nodes in the 
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To see what 



tree would be the sum of for i from 1 to 8, 
the upper and lower bounds on are# consider again the 
tree in fig. 11. The fewest number of children of each node 
is found in the leftmost path down the tree. In this pathr 
every node has two children. In the rightmost oath down the 
tree/ a node at level i has i children and this is the larg- 
est number of children for any node at that level. Thus the 
total number of nodes at any level i lies somewhere between 
the two cases/ or 

2 " <= 0^ <= n ! 

Since 8^ is growing faster than 2*^ it is easy to see 
that as n gets large/ the possible size of the tree is going 
to far outweigh any of the polynomials so far discussed. 
Obviously/ unless a significant amount of pruning is done/ 
the time to expand every node at each increasing level in 
the tree is going to grow far too rapidly. Since the search 
strategy was fixed for all orderings/ the most significant 
way to compare the time savings of each ordering was to 
count the number of nodes which needed to be expanded. 
Therefore/ the number of calls to the function EXPAND was 
used as the primary measure of the time complexity of the 
various orderings. 
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2. Space 



There are only two data structures in these programs 
whose size is dependent upon the vertex ordering: the array 
of nodes in the search tree and the array of pointers used 
for the priority queue . But the number of nodes in the 
queue is always less than the total number of nodes in the 
tree since expanded nodes are removed from the queue. Since 
the number of expansions has been recorded for a time 
analysis. only the total number of nodes in the search tree 
at the time an optimal solution is found need be recorded 
for a space analysis. 

In order to see the shape of the search tree, anoth- 
er parameter, the average branching factor, was recorded for 
each level in the tree. The average branching factor is de- 
fined to be the average number of children of each node in 
the search tree. A full binary tree would thus have an 
average branching factor of 2.0. If the branching factor is 
recorded for each level in the search tree, one can see the 
'shape* or the 'narrowness' of the tree. The most desirable 
situation is the one in which the branching factor stays at. 
or as close as possible, to 1.0 in the first several levels 
of the tree. This indicates that the maximum pruning took 
place as early in the search as possible, thus reducing the 
search space by more nodes for each prune. The 'goodness' 
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of the ordering functions can thus be seen by how long the 
branching factor stays at 1.0 and how low it stays when it 
does start to increase. 

B. TIME ANALYSIS 

In analysing the time required by the program for each 
ordering method there were two cases of interest: average 
timer and the time required in the worst case. Figure 12 is 
a plot of the average number of expansions required for each 
ordering method when graphs containing 30% of the maximum 
edges were colored. 

For the worst case graphs the natural logarithm of the 
number of expansions was plotted since exponential growth 
was expected. Figure 13 is a plot of the worst case graphs 
for 30% edges. These plots are not smooth and in fact are 
quite ragged because each point actually represents the 
coloring of one particular graph in the sample of 100. Note 
that the ordering methods parallel each other in the dips 
and hills of the plot/ indicating that the same graph was 
usually the most difficult to color for all of the methods. 

The significance of these two plots is in the total dom- 
inance of the LEAST COLOR AVAILABLE method of ordering the 
vertices. Though all the other methods eventually 'explod- 
ed' even in the average case/ LCA stayed very close to 
linear for the full range of graphs up to 32 vertices. 
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Figure 12 
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LN (MAXIMUM EXPANSIONS) 




MAXIMUM EXPANSIONS VS. GRAPH SIZE 
Figure 13 
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Fur t he rmo r e # the worst case analysis shows the number of ex- 
pansions for LCA leveling off between 20 and 30 vertices. 
It would be an interesting experiment to carry the graph 
sizes out past 32 to see if and when LCA also starts to show 
explosive growth. 

Another result/ though not shown in graphical form is 
that the time requi red by any of the ordering methods is 
lowest for graphs with very few edges as well as for graphs 
which are close to being complete. For any given value of 
n, empirical evidence indicated that the most time is re- 
quired when the number of edges is in the vicinity of 50% 
(see appendix B). This is an intuitively obvious result if 
one considers what happens in the search tree. With very 
few edges/ very few colors will be needed so the number of 
children of a node will be minimal. The search/ in this 
case/ will tend to proceed down the left side of the search 
tree. Most of the right side will be pruned away because 
the cost of the nodes in the right side of the tree will be 
greater than the cost of the optimal solution. 

In the case where there are many edges/ a new color will 
be required at each level and the search will tend to stay 
very close to the right side of the tree. The left side 
will be pruned away since many of the nodes on that side 
will represent illegal colorings due to edges. 
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C. SPACE ANALYSIS 



Just as for time/ the space analysis is of both average 
and worst case situations. Figure 14 is a plot of the aver- 
age number of nodes in the search tree when the first op- 
timal solution was found. Figure 15 is a plot of the natur- 
al log of the maximum nodes needed for any of the 100 
graphs. Again/ the number of edges is 30% of the maximum 
poss i b 1 e . 

Not su rpr i s i ngl V / the space analysis results are very 
similar to the results found in the time analysis. In every 
case LCA is superior in the sense that it uses less soace 
and its growth rate is much improved over the other methods. 

Since the relationship between the number of expansions 
and the total number of nodes is the branching factor/ it 
comes as no surprise that LCA exhibits the lowest branching 
factor at the top of the search tree. Figure 16 shows the 
average branching factor as a function of depth for each of 
the 4 heuristic orderings. Note that LCA and ADJACENCY have 
the same branching factor for 3-4 levels but then ADJACENCY 
grows faster in the midrange. Note also that both of these 
methods have a higher branching factor in the middle of the 
tree than either DEGREE or DVEC 3. This fact shows how 
oowerful 3 or 4 prunings at the top of the tree can be. 
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AVERAGE BRANCHING FACTOR 




AVERAGE BRANCHING FACTOR VS. DEPTH IN TREE 

Figure 16 
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D. ADVANTAGE OF DIVIDING THE GRAPHS 



In the first stage of this coloring process# the graphs 
were colored just as generated with no attempt to divide 
them into biconnected components. In the second stage# they 
were divided and then each component was colored and the 
number of expansions reduired was summed up for all the com* 
ponents. For sparse graphs where there were many very small 
components# the results of this experiment were quite sub* 
stantial. However# as the graphs became more dense# the 
original graph was biconnected and no savings were derived 
f rom dividing. 

Figure 17 shows the sire of the largest biconnected com* 
ponent as a function of the number of edges. In this plot# 
the ratio of the size of the largest component to the origi- 
nal graph size is plotted for 4 cases. Note that since the 
maximum number of edges grows as n^ , when the number of 
vertices in the graph increases# the graph becomes bicon- 
nected with a smaller percentage of edges. Thus it is rea- 
sonable that the empirical evidence gained here indicates 
that a graph with 100 vertices is biconnected 99% of the 
time with 6% edges whereas a graph with 50 vertices is 
biconnected 99% of the time when it has 12% of the maximum 
vertices included. 
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In investigating the effects of dividing on graph color- 
inpf note from fig. 17 that with 30% edges# the largest 
biconnected component is very close to the original graph 
for all n. Therefore# in order for dividing to be benefi- 
cial# random graphs must be very sparse. For the empirical 
results shown in fig. 18 graphs were generated with from 10% 
to 30% edges. The benefit of dividing is then shown by 
plotting the ratio of the number of expansions needed for 
the divided graph to the number of expansions for the undi- 
vided graph against the percentage of edges present. 
Inasmuch as the relative benefit was the same for all order- 
ing methods# only LCA was used for this graph. 

It appears as though the benefits from dividing are in- 
significant since only very sparse random graphs have arti- 
culation points. But many real-life problems which are 
modeled by graphs may be quite dense yet have articulation 
points. This is true particularly in transportation and 
communications problems where some of the vertices may 
represent junctions or hubs through which all traffic must 
flow. Consider the city of Chicago in a nationwide tran- 
sportation network for example. In these cases graph divid- 
ing may result in a much faster coloring. 

Also there are problems which have graph models which 
are indeed very sparse. Consider# for example# the exam 
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scheduling problem in a large university. The vertices of 
the grapn are the classes and there is a edge between two 
vertices if at least one student is enrolled in both 
classes. The question then is: How many exam periods are 
needed to schedule all final examinations with no conflicts? 
Consider the hypothetical but realistic case of a university 
with 200 classes/ and a average of 20 students per class. 
Since many students take the same curriculum together it is 
not inconceivable to have a graph with an average degree of 
each vertex of 5. This would result in a graph of 200 ver- 
tices and 500 edges or 2 1/2 %. This is a very sparse graph 
which is likely to have several articulation points. In 
fact/ the biconnected components are likely to be the vari- 
ous departments of the university. Dividing the graph prior 
to solving this problem would surely result in a time and 
space savings. 
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IV. CONCLUSIONS 



The branch and bound algorithm described in this oaper 
can be used to find the chromatic number of graphs. Though 
no empi rical evidence was gained for graphs larger that 32 
vertices# it is felt that somewhat larger graphs can also be 
colored in a reasonable amount of time on the average with 
the LEAST COLOR AVAILABLE method of ordering the vertices. 
Further# some improvement in the lower bound function or 
search strategy might be possible which would make the algo- 
rithm more efficient. For example# some lookahead procedure 
might conceivably improve the method of assigning costs to 
each node of the search tree. Also# the priority of a node 
in the gueue was primarily the number of colors in the par- 
tial coloring and only secondarily a function of the depth 
in the tree. A possible improvement is to make the priority 
some algebraic combination of the two factors thus giving 
more emphasis on the depth in the tree. This might result 
in some faster colorings. Since the emphasis here was on 
finding the best branching function# little effort was ex- 
pended in looking for improvement in the other 2 components 
of the branch and bound strategy. 

For the five branching functions investigated# LCA seems 
far superior for all graphs in both the time reouired and 
the amount of computer storage space needed. Though all 
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five functions have worst case time and space complexity 



which appears to be exponential/ the LCA method has an aver- 
age case time complexity which appears to be nearly linear 
in the range of graphs from 5 to 32 vertices. 

Though the research and experiments done here focused on 
the exact algorithm for graph coloring, a few very minor 
changes would result in an approximation algorithm. The 
function EXPAND could be changed such that it only generates 
the first possible child node rather than all children of a 
node. No priority queue would be needed and the number of 
nodes required in the search tree would just be n, the 
number of vertices in the graph. So the space-complexity of 
the algorithm would be linear in n. The time required would 

4 

be 0(n ) if LCA were used for an ordering function/ but it 
would be interesting to see how good its ability to approxi- 
mate the real chromatic number is compared to other order- 
i ngs . 

Dividing a graph at its articulation points prior to 
coloring is an effective way to reduce the time required to 
find the chromatic number. However/ unless the graph is 
known to have articulation points or unless it is very 
sparse/ the effort spent in finding the biconnect ed com- 
ponents may not be worth it. When a graph does have articu- 
lation points/ there is a two-fold benefit gained by div id- 
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i ng the graoh. The first is/ of course^ that 
the resulting graph is smaller and thus easier 
The second benefit comes from the fact that the 
components tend to be dense which makes coloring 



the size of 
to color, 
bi connected 
eas i er , 
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APPENDIX A 



This appendix contains the source listing for the com- 
puter programs used in this research. The program is organ- 
ized into various functions which share some common global 
variables. These variables are listed in the module EXTERN 
DECS. The remainder of the modules contain one or more 
functions grouped together according to their logical func- 
tion, The modules and the function in this appendix are: 



MODULE 

EXTERN DECS 
MAIN 

GENGRAPH 

ORDERNODES 

COLORALL 

DIVIDE 

TREEMANAGER 



FUNCTIONS 

none 
ma i n ( ) 
gengra f ( ) 
random ( ) 
ordernode ( ) 
p i c knode ( ) 
CO 1 ora 1 1 ( ) 
expand ( ) 
convb ( ) 
d i V i de ( ) 
proc ( ) 
getnode ( ) 
addq( ) 
get t op ( ) 
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/»***************************** 

GLOBAL MACRO DEFINITIONS 



n * * * 


***** 


* * * 


★ ★ 


ififitlfUififititif 


if it it if it if / 


«de f i ne 


NUMNODES 


aooo 


/★ 


Maximum number of 


tree nodes */ 


»de f i ne 


NUMQ 


2500 


/ ★ 


Maximum size of the priority Q */ 


i*def i ne 


NULL 


0 


/★ 


The empty pointer 




f i ne 


NUMTRIALS 


100 


/if 


The number of trials for each size 










of graph */ 




#def i ne 


MAXSIZE 


30 


/★ 


The largest graph 


size (n) */ 


#de f i ne 


MINSI2E 


10 


/if 


The smallest graph 


*/ 


/ * * * 


****** 


* * * 


★ ★ 


★ ★★★★★★★★ 


it it H if if it if 



EXTERNAL DECLARATIONS 

********<>**j>****i>**<>****i***** * / 



/* trnode is the structure of the nodes 
which make up the search tree */ 

struct trnode < 

unsigned depth :8; 
unsigned cost :0; 
unsigned vertex :8; 
unsigned color :0? 
struct trnode ^parptr; 

>; 



/* tnode is a array of search tree nodes, 
q is an array of oointers to tree nodes */ 
struct trnode t node (NUMNODESj ^ *q INUMQJ ; 

/* ndata is a record of data regarding the 
graph vertices to be used in ordering the 
vertices by the branching function, 
order is an array of records */ 



struct ndata { 

unsigned next 
unsigned degree 
unsigned flag 
> order 132J } 
i nt f d? 
i nt n; 
int count/ 
i nt hi node / 
int index; 



5; 



5; 

i; 



int 
i nt 
i nt 
i nt 
int 



qs 1 ze; 
c h rom ; 
color C32) 
ordnum; 
saven ; 



/* file descriptor for output */ 

/* number of vertices */ 

/* number of calls to EXPAND */ 

/* first vertex to color, */ 

/* index of next available node in the 
array of tree nodes • also the number 
of nodes in the search tree at any time 
index into the priority queue */ 

The chromatic number of the graph */ 
The color class of each vertex */ 



/* 

/* 

/* 

/* 

/* 



The ordering function being used */ 
The saved value of n for DIVIDE */ 



*/ 
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int gra f (321 132] ; /* The adjacency matrix */ 

float numsubs; /* The number of biconnected subgraphs */ 
float bigsubn; /* The size of the largest biconn, como, */ 
float bigsube /* Number of edges in bigsubn */ 

long int graph C323 ; /* The adjacency matrix A in bit strings */ 

int sonct 1321 ; /* A counter of child nodes for each level */ 

int parct(32i; /* Parent counter for each level */ 



/****************************** 

MAIN 

***************************** */ 

/* This function provides for the flow of 
control throughout the program and records 
the data which is to be collected. It also 
contains the code for any output generated */ 

main ( ) { 



i nt f ul 1 » /* 

int edges; 
i n t i » j ; / * 

int ndf 1 ag [6] ; /* 

int trial; / * 



float maxedges; /* 
float maxbrCbJ; /* 
float maxndtbl; /* 
float avgbr £6] ; / * 
f 1 oat avgnd £61 ; / * 
char * 1 abel £61 ; 



The percentage of edges */ 

The number of edges */ 
iteration variables */ 

Too large of a search tree flag for each 
ordering method */ 

The trial number for each graph size */ 
The maximum edges possible for graoh */ 
The maximum expansions required */ 

The max nodes required */ 

The average expansions required */ 

The average nodes in the search tree */ 



The outDut file and labels 
fd = copen ( " /wor k/reraut en/out 3" » ' a ’ ) ; 



1 abe 1 
1 abe 1 
1 abel 
1 abel 
1 abel 



£11 

£21 

£31 

£ai 

£51 



"RANDOM"; 
"DEGREE"; 
"DEG 3"; 
"AOJCEN"; 
"LCA"; 



*/ 



/*. For each value of n desired */ 
for(n=MINSIZE;n<=MAXSIZE;n+=5) { 

/* compute maximum edges possible */ 
maxedges = (n * (n-l))/2; 

/* For each value of the percentage of 
edges des i red */ 
for(ful 1=2; fulK = 20; ful 1 t = 2) { 
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/* comoute number of edges */ 
edges = full * maxedges/ 1 00 ; 

/* Print headings */ 
orintf(fd/"0 =X3d» EDGES =%4d ( %ad%% ) 0 , 
n , edges / full); 

printf(fd/"%5s%10sX10s0»"order"/"branches’'/"nodes"); 
pri nt f ( fd/ "X15s212s" # "avg max"#'*avg max "); 
print f(fd#"0); 

/* Initialize data collection variables */ 
for(i=l;i<=5;i++) { 
maxb r I i ) = 0 ; 

maxnd ( i ] = 0 ; 

avgbr til = 0 ; 
a vgnd til = 0 ; 
ndf lag ti 1 = O; 

> 

/* For the number of trials desired */ 
for(tria1=i;trial<=NUMTRIALS;trial++) < 

/* Generate a random graph */ 
gengraftn/edges); 

/* For each ordering method */ 
for(ordnum = 5;ordnum< = 5;ordnum + 'f) { 

if(l ndf 1 ag lordnum] ) { 

/* Divide the graph and color the components * / 
if(i dividet)) ndf 1 ag tordnumJ = l; 

/* Record data */ 

avgbr tordnuml += count; 
a vgnd to rdnuml += index; 
iftcount > maxbr tordnuml ) 
maxbr tordnuml = count; 
iftindex > maxnd tordnumJ ) 
maxnd tordnuml = index; 

> 

> 

> 

/* Print out the data collected */ 
f or ( i = 1 ; i <=5 ; i ++ ) { 

i f ( i ndf lagtil ) { 

/* Compute the averages */ 
avgbrtiJ /= NUMTRIALS; 
avgndtil /= NUMTRIALS; 
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> 



> 

cex i t ( ) ; 



} 



printf(fd/'*%6s"»1abel [i]); 

printf(fd/"Xy.0f%5.0f",avgbrCi] rmaxbrtil ); 
printf(fd#"%6,0f%5,0f "favgndtiJ »maxndli] )J 
printf(fd»"0); 

> 

else printf (fdr"X6s NOT ENOUGH SPACE AVAILO, 
label [ n ) ; 



1 



/****************************** 

GENERATE GRAPH 

**********★*******************/ 
/* This function generates a random graph with 
n vertices and e edges and records it in a 
global adjacency matrix */ 

gengraf(n»e) int n#e? < 

i n t i > j r a f b ; 
long int convb()» 

/* initialize the adjacency matrix */ 
for ( i =0 ; i < = n ; i ♦ + ) { 

graph ( i J = 0 ; 

f or ( j =0 » j <=n ; j ++ ) graf til Cjl = 0? 
grafCiltil = 1; 

} 

while (e > 0) { 

/* Get 2 random vertices */ 
a = random(n); 
b = random (n ) ; 

/* Find the next ungenerated edge */ 
while(graf[aJ[bl == 1) { 

a=(a%n)+i; 
i f (a == 1 ) b = (b 2 n) + 1 ; 

} 

/* Record the new edge */ 
graf Cal Cbl s 1 ; 
grafCblCa) = 1; 

e--; 

> 

/* Put adj matrix in array of bit strings */ 
for (i=l;i<=n;i++) { 



aa 



I 



I 

I 



graft nCi) = 0; 

for(j=i;j'^ = nfj'*"*’) iftgraftnCj]) graphti] '= convb(j); 

> 



/t ***************************** 

RANDOM 

***************************** */ 
/* This function is used to get random numbers 
between 1 and the input argument n. 
taken from Grogono [ref. *”J */ 

int random(n) int nt { 

static long int seed 18403; 
inti; 

seed = (25173 * seed + 13849) % 65536; 
i = seed X ni 
return( i t 1 ) ; 

> 



/*****************it************ 

OROERNODE 

***************************** */ 

/* This function orders the vertices of 
the graph for each of the static ordering 
methods. It does so by filling the array 
order. next with the vertex to follow each 
vertex in the graph. */ 

ordernodeC) { 



int i / i » Ic ; 
int save; 
int t emph i ; 
int max; 
int col sum C32) 



i nt 
i nt 
f 1 oat 
f 1 oat 
1 ong 
1 ong 



/* iteration variables */ 

/* the saved last vertex */ 

/* the temporary next vertex */ 

/* any maximum needed to maintain */ 

/* used in Adjacency ordering to keep 
track of number of previously colored 
vertices a vertex is adjacent to. */ 

The column with the max I's in colsum 
The highest degree of any vertex 
/* The vector of degrees in DVEC 3 */ 
The max vector calculated in DVEC 3 */ 
mask; /* A bit manipulation mask */ 
convb(); / * A converting function (SUPP MOD) * / 



maxcol ; /* 

hideg; /* 

dvec (4) (32) ; 
maxdvec; /* 
i nt 
i nt 



*/ 



max = 0; 
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f or ( i 



/* initialize the vertex data * / 
0;i<=nJit+) < 

order [ i ]. degree = 0; 
order(i].next = 0 ; 
order C i J . f 1 ag = 0 ; 
mask = graph til; 

/* Compute the degree of each vertex */ 
for(j=l?j<=n?)++) { 

if ((mask & 01) == 1) o rder [ i J , degree += l; 
mask >>s 1/ 

> 



/* find the max degree and the vertex */ 
i f (order (i ) .degree > max) { 

max = orde r ( i 1 . degree ? 
h i node = i > 

> 

> 



/* Order the vertices depending on the 
current ordering number */ 
swi tch (ordnum) ( 
case 0 : 

break ; 

/* RANDOM */ 

case 1 : 

f or ( i =0 ; i <n; i + + ) order (i ] .next = i t 1; 

hinode = i; 

break; 

/* by DEGREE */ 

case 2: 

order (OJ .next = hinode; 

/* .flag indicates already colored vertex */ 
order (h i node) . f 1 ag = i; 
save = hinode; 

/* For each vertex find the next one */ 
for(i=l;i<=n;it+) { 
max = 0; 



/* Find the highest 
vertices */ 
for(J=i;j<=n;jt+) { 

if((orderlj] .degree 



degree of the uncolored 



> max)&&(! order ( j ]. f 1 ag) ) 



{ 
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max = order [ n .degree; 
temph i = j ; 

> 

> 

record the next vertex */ 
order Isave) .next = temphi; 
order ( t emph 1 J . f 1 ag = 1; 
save = temph i i 
> 

break ; 

/* by OVEC 3 */ 

case 3; 

/* initialize 0th iteration to degree */ 
f or ( i =1 ; i < = n; i ♦+ ) dvec(0](i] = order ti ] .degree; 

/* For 3 iterations */ 
f o r ( i = I ; i <=3 ; i ♦+ ) { 

/* For each vertex */ 
f or ( J = I ; j < = n ; j -tt ) { 
dvec ( i ] t i ] = 0 ; 

/* For every other vertex in the graph */ 
f or ( k = 1 ; k<=n ; k t + ) { 

mask = convb ( k ) ; _ 

/* If they are connected compute summation 
for the ith iteration */ 

if((graph[j] i mask) > 0) 

dvec(i) Ij] += dvec(i-l) IkJ ; 

> 

> 

> 

/* Order the vertices */ 

/* The ordering method is identical to case 2 
except for the ordering parameter */ 
save = 0; 

for(i=l;i<=n;i++) { 

maxdvec = O; 
for ( j =1 ; j <=n; j + + ) { 

i f ( (dvec (3) n 3 > maxdvec) && (i order I j 3 . f 1 ag ) ) ( 
maxdvec = dvec 133 Cj3; 
t emph i = j ; 

> 

} 
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order (save] . nex t = temphi; 
order ( t emph i ]. f 1 ag = 1; 
save = t emph i ; 

> 

hinode = order (0) , nex t ? 
break ; 



/* by ADJACENCY */ 

case 4: 

/x initialize */ 
f or ( i =l ; i <2n; i ♦+) co1sum(i] 
order (01 .next = hinode; 
order (h i node! . f 1 ag = 1; 
save = hinode; 
colsumisavel 33; 



o; 



/* For each vertex find the next one */ 
f o r ( j = 1 ; ] <=n- 1 ; i ♦ ♦ ) { 
max = O; 
mask = 01; 

/* increment the bit string colsum with each 
vertex added to the set of colored vertices */ 
f or ( i = I ; i < = n ; i +■♦• ) ( 

i f ( (graph (save! & mask) > 0) 
col sum ( i 1 ; 

if(colsum(i] > max)max = colsum(iJ; 
mask <<= i; 

} 



/* Order the vertices - again the same 
ordering method is used as in case 2 except 
the ordering parameter is by value of colsum */ 
hideg = O; 
f or ( i = 1 ; i <=n; i t+ ) i 

i f ( (col sum (i ) = = max)8t&(order (i ] .degree>hi deg) 

&&(! order ( i ). f 1 ag ) ) ( 
hideg = order (i ] .degree; 
temph i = i ; 

> 

> 

order (save] .next = temphi; 

order (temph i ]. f 1 ag = i; 

col sum (temoh i ] -= 33; 

save = t emph i ; 

> 

break; 
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/* by LEAST COLORS AVAILABLE 

case 5: 

no work done in this function. 





See 


PICKNODE */ 








> 

> 


break ; 










/ * * * 


★ ★★★★★★ 
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k k k 


★ ★ 
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PICKNODE 








* * * 


it if if -k it -k 




k k k 


k k 


k k k k k k 




This 


f unc t i on is called 


from 


the 


rout i me 




EXPAND to get the next 


vertex 


t 0 


color. 



If one of the static orderings is being 
used it just looks uo the next vertex in 
the table order.next. If the dynamic 
ordering LCA is being used/ it must trace 
the search tree to the root to determine 
the next vertex to color, */ 

pi c knode ( nodept r ) struct trnode ^nodeptrJ { 

/* Declarations */ 
long int cc [321 / vdone/ 1 ongl ; 
int maxfisum/k/i,sumf32) /nodenum/hideg; 
struct trnode *tempptr/ 

/* If static ordering/ just table lookup 
if(ordnum < 5) ret urn ( order Cnodept r->vert ex J , next ) / 

/ * initialize 

vdone is a bit string to record which 

vertices have already been colored, 

cc is the color class matrix 

sum counts the colors not avail for each 

vertex, 

vdone = O; 
f or ( i = 1 / i < = n ; i -f-t) { 

cc li 1 = 0; 
sum I i 1 =0 ; 

> 

temoptr = nodeotr? 

/* Trace from node N to root, */ 
while (tempptr 1= NULL) < 

fill color class matrix */ 
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cc [temppt p->col or] 1= graph [temppt r->vertex] ; 

/* record colored vertices */ 
vdone 1= convbCtemopt r->vertex ) » 
tempptr = temppt r->parpt r» 

} 

1 ong 1 = 01; 

max = O; 



/* For each of the n vertices */ 
for(k=i;k<=n;kt+) { 



/* If the vertex k is not colored */ 
if((vdone & convb(k)) 0) { 

i sum = 0 ; 



/* count the color classes to which vertex 
k cannot be added */ 
f or ( i = I ; i <snodept r->cost;i { 

i sum +■= longl i ( cc C i J >> ( k - 1 ) ) ; 

> 

sum (kJ = i sum; 

/* keep track of fewest number of color 
classes available */ 
i f ( i sum > max) max = isum; 

> 

> 

hi deg = 0; 

/* Of all the vertices with the fewest color 
classes available^ find the one with the 
highest degree */ 
f or ( i =1 ; i < = n; i ■>■•♦■) { 

if((sumCi} == max) && (l(convb(i) & vdone))) { 
i f (order { i J .degree > hideg) ( 
hideg = o rder I i }. degree ; 
nodenum = i ; 

} 

> 

} 

/* Return the vertex number of the next 
vertex to color. */ 

return(nodenum) ; 

} 
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/★ ★ A ★ A 



COLOR ALL 

■k'kif'ki^itiK’k-kififiX'k'k'kif’kifiic-kiiitiKi^'kiKiKififit 

This function creates the root node of the 
search tree/ puts it in the priority queue 
then calls the function EXPAND to expand 
the top node in the queue until the top of 
the queue is a complete coloring, */ 

coloralH) { 

struct trnode *nodept r / *get top ( ) ; 
i nt i ; 

/* Initialize for new coloring */ 

k = o; 

index = 0; 
qsize = O; 
count = 0; 
tnode tO] .cost = 31; 

for ( i =0; i <NUMQ; i ++) qtil = StnodelO); 



/* Get an empty node/ fill it with the 
data for the root of the tree and 
out it in the queue, */ 
if (i (nodeptr = getnodel))) return(O); 
nodept r->parpt r = NULL; 
nodept r->dept h = i; 
nodept r->cost = l; 
nodept r->vert ex = hinode; 
nodept r->co 1 or s i; 
addq(nodept r ) ; 

/ * 'forever' * / 



f or ( ; ; ) { 
nodept r 



f* Get the too node in the queue */ 
get top ( ) ; 



/* If this is complete coloring/ quit */ 
i f ( nodept r->dept h == n) break; 
else if (1 expand (nodept r ) ) return(O); 

> 



/* Record the chromatic number */ 
chrom = nodept r->cost ; 

/* Record the color class assignment of 
each vertex in this oath, */ 
while (nodept r->parpt r i= NULL) ( 

col or (nodept r->vertex] = nodept r->col or ; 
nodeptr = nodept r->oarpt r ; 
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I 




nodept r->co lor; 



> 

C0 1 OP Inodept r»>vert ex 3 = 

retupn(l); 

} 



/*********★******************** 

EXPAND 

*J>***T************************* 

This function expands the node pointed to by 
the variable nodeptr. It creates a child 
node for each possible color available as 
well as a node for one new color, */ 

expand (nodeptr) struct trnode *nodeptr; { 

struct trnode *temppt r, ^newnoder *getnode ( ) ; 
i nt vert exnum^ i ; 

long int mask ^ co 1 ava i 1 / convb ( ) ; 

/* increment counters */ 

count ♦= i; 

parct (nodept r->dept hi ♦+; 

/* Get the number of the next vertex to color */ 
vertexnum = p i c knode ( nodept r ) ; 

/* create new nodes for all colorings using 
only existing colors */ 

CO 1 ava i 1 • 01; 

for ( i = 1 ; i <nodept r->cost ; i ++ ) < 

col avai 1 <<= 1 ; 
colava il }= 01; 

/* colavail is a bit string with a one in the 
position of every available color, */ 

> 

tempptr = nodeptr; 

/* For each node N in the path to the root */ 

for ( ; ; ) ( 



/* If N, vertex and vertexnum are connected */ 
i f ( (graph (vertexnuml & convb ( temppt r->vertex ) ) 1= 0) { 

mask 5 convb ( t emppt r->co 1 or ) ; 



/* Remove N, color from colavail */ 
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if ((colavail & mask) 1= 0) colavail t= mask; 

> 

if ( temopt r«>parot r == NULL) break; 
tempptr = temppt r«>parpt r; 

) 

/* If any colors are left */ 
i f (col avai 1 1=0) { 

mask = 01; 

/* For each color available */ 
for ( i =l ; i <=nodept r->cost ; i ++) < 

if ((colavail & mask) != 0) ( 

/* Generate child node */ 
if ((newnode = getnode()) == 0) return(O); 
newnode">dept h = nodept r->deot h + i; 
newnode*>cost = nodept r->cos t ; 
newnode">ver t ex = vertexnum; 
newnode->co 1 or = i; 
newnoOe->oarpt r = nodeotr; 

/* increment child counter */ 
sonct (nodept r->dept hi + + ; 

/* Put the node in the queue */ 
addq(newnode) ; 

> 

mask =<< i; 

> 

> 

/* Generate one more node with a new color */ 
if ((newnode = getnode()) == 0) return(O); 
newnoae*>dept h = nodept r->dept h + l; 
newnode">co5t = nodept r->cos t t i; 
newnode~>ver t ex = vertexnum; 
newnode->co 1 or = nodept r->cost + l; 
newnode*>parpt r = nodeptr; 
sonct (nodept r->dept hi +■^; 
addq(newnode) ; 
return ( 1 ) ; 

> 



/*H*it1tit1titit***1t**it*****if*it-M***itit 

CONVERT TO BIT 

**★*******♦**★****************/ 

/* This function converts an integer argument 
between 0 and 16 to a bit string which has a 
1 in the position of the argument and O's 
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everywhere else */ 
long int convb(arg) int arg» { 

i nt i ; 

1 ong int bi tpos ? 
bi tpos = 01; 

for(i=l/i<arg;it+) bitpos =<< i; 
return(bi tpos) ; 

} 



DIVIDE 

/* This function divides the graph into its 
biconnected components. It then creates a 
new adjacency matrix far each component and 
calls the function COLORALL to color the 
component. It records the data for each 
component and sums it up for all comps of a 
graph. When it returns/ the sum data is in 
the variables set up by MAIN, */ 



divide!) { 



i nt 


numbe r [321 ; 




i nt 


back 132 ) f 


/★ 


1 nt 


St ac k V (32) i 


/* 


1 nt 


num ; 


/it 


i nt 


w / V » t op ; 


/it 


1 nt 


St vpt r ; 


/ it 


i nt 


St ept r ; 


/it 


int 


nume i 


/it 


i nt 


numv ; 


/it 


i n t 


vnum [321 i 


/it 

t 


1 nt 


totcount; 


/* 


i nt 


tot i ndex ; 


/it 


i nt 


table (32) ; 


/it 


i nt 


t ab 1 ept r ; 


/it 


1 nt 


vert 1 / vert2; 


/it 



long int 
int i / j ; 



a back 



edge exist * / 



*/ 



components */ 



convb ( ) ; 



/* An edge in a biconnected component defined 
by the two vertices at its end */ 
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struct edge { 

uns i gned v 1 : 8 > 
unsigned v2: 8? 

> s t ac ke [3501 ; 

/* Stacke is an stack of edges 

/* Save off the oringginal value of n - the 
biconnected components will be different */ 

saven = n; 

/» Initialize all variables */ 

number 1 1 ] = 1 ; 

num 2 2f 
stackvtOl s 1; 
too = I ? 

f OP ( v=2 ? v<-n; v+ t ) numberlvl = 0; 

steptr = -i; 
stvptr = 0; 
totcount = 0; 
totindex = 0; 

/* 8egin a depth first search of the graph 
exploring forward edges */ 

forward : 

/* While there is an unexplored forward edge 
from the top vertex on the stack to some 
vertex w */ 

while((w = proc(top)) 1= 0) { 

/* Stack the edge (top/w) on the edge stack */ 
s t ac ke I + + st ept r] . V 1 = top; 
stacke (steptrl ,v2 = w; 

/* If w has already been explored give the min 
of number(w) and back(top) to back(top)* */ 
i f (number twl > 0 ) { 

if(number(wJ < backCtopl) backttopl = number (wJ ; 

> 

/* Else stack w and explore it */ 

else { 

number [wl = num++; 
backtwl = numoer twl ; 
stackv [f + st vpt rl = w; 
top = w; 

> 

> 



/* If there is no forward edge then backup 
to the last vertex which still has a forward 
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edge which is unexplored */ 



bac kup : 

i f (st vpt r >0) { 

/* try the next to top vertex */ 

V = St ac kv [st vpt r - IJ; 

/* if it has a number less that the back of 
the too vertex then we have an articulation 
point and the bi connected component is on 
top of the edge stack */ 
if(backttop] >= numberfvl) < 

/* increment and initialize counters * f 
numsubs l; 
nume = 0; 
numv = 0; 

for(i=l;i<=saven;it+) graphti] s 0; 
tableptr = 0; 
tablelO] = O; 

f or ( i =0 ; i <=sa ven ; i +t ) vnumtil = 0; 

/* POP edges until end of component */ 
wh i 1 e ( ( numbe r I s t ac ke [s t ept r 1 , V I ] >= numberfvl) 

(number [stacke [steptr] .v2J >= numoer [vj ) 

(steptr >= 0)) { 

/* increment data counters */ 
nume + ♦ ? 

i f ( + + vnum (st ac ke (s t eot r 1 . V I ] == 1) numv + + ; 

i f ( + + vnum (s t ac ke (s t eo t r J . v2] == 1) numv + + ; 

vertl = 0 ; 
vert2 = 0 ; 

/* check to see if either vertex is already 
in the table - if so» get the index which 
is going to be the new vertex number */ 
f or ( i =0 ; i <=t ab 1 ept r ; i ++ ) { 

if(table[il == st ac ke ( st ept rl , v 1 ) vertl = i; 
if(table(i] == stactce (stept r] . v2) vert2 = i; 
} 

/* if not in tables put it there and get new 
vertex number 
i f ( ve rt 1 == 0 ) { 
t ab 1 ept r + + ; 

table Itabl ept r] = st ac ke (s t ept r J , v 1 ; 
vertl = tableptr; 

> 

i f ( vert2 == 0 ) { 

t ab 1 ept r + + ; 

table (tableptr] = stacke (stept rJ .v2; 
vert2 = tableptr; 

} 
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/* record the edge between the new vertex 
numbers */ 

graphCvertll 1= convb ( vert 2) ; 
graph {vert2] 1= convD ( vert 1 ) / 

/* decrement edge stack pointer and do again */ 
stept r--; 

> 

/* The biconnected component is now recorded 
in the adjacency matrix and ready to color */ 

/* Is this the biggest component of this graph */ 
ifCnumv > bigsubn) bigsubn = numvj 
if(nume > bigsube) bigsube = nume; 

/* new n is size of the table */ 
n = tableptr? 

/* Call the ordering function */ 
ordernode ( ) ; 

/* Color the comoonent */ 
if(i coloralH)) return(O); 

/* Increment the cummulative counters */ 
totcount += count; 
totindex t= index; 

> 

/* If no art. pt. pop the vertex stack and 
explore any forward edges of the vertex 
which is now on top */ 

else if(back[top] < backtvl) backCvl = backltopj; 
top = St ac kv (--St vpt r J ; 
goto forward; 

} 



/* When all the connected vertices have 
been explored# check for another unconnected 
comoonent and explore it */ 
f or ( i =2 ; i <=saven ; i t+ ) ( 

i f (number ( i J == 0 ) < 

number ( i 1 = 1 ; 

num = 2; 
stackviOl s i; 
top = i; 
stvptr = 0; 
stept r = - 1 ; 
goto forward; 

} 

} 

/* When everything is done and all components 
are found and colored# set the counters to the 
cummulative totals# reset n# and reset the 
adjacency matrix */ 
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i 








count = totcount/ 
index = totindex; 
n = seven? 

for (i=l/i<=n?i++) { 

graf [i ) t i J = 0 ; 
for (j=l;j<=n;jt+) { 

if (grafliJljl == 2) graftiJtjJ = I? 

if (graf {ill]] == 1) graph(i) 1= convb(j); 

> 

> 

return(l); 

} 



/**************************<>*** 

PROC 

*****♦**♦*******♦*************/ 

/* This function determines if we can 
proceed in the depth first search exploring 
forward edges. It marks the adjacency 
matrix with a 2 as each edge is explored 
so that no edge is explorea twice */ 

proc ( i ) i nt i ? ( 
int j; 

/* Is there a vertex adjacent to vertex i 
on an edge which hasn’t been explored. If so 
mark it explored and return the vertex number 
for( j=i; j<=saven; j++) 

i f (graf lil Ijl == 1) { 

graf (i 1 Ij) = 2? 
graf {jl li) =2; 
return(j); 

> 

return ( 0 ) ; 

} 



/****************************** 

TREE MANAGER 

itHit-k**** -ft Itititliitit***************/ 

/* This module contains the functions 
necessary to manage the search tree nodes 
as well as the priority Queue. */ 
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GETNOOE 

KK****************************/ 

/* This functions maintains the array of 
nodes and returns a pointer to the next 
available node when called. Index is a 
global variable */ 
struct trnode *getnode() { 

/* index gets initialized in COLORALL */ 

index += 1? 

if (index < NUMNODES) return(&tnode (i ndex] ) ; 
el se return(O) ; 

) 



/*************************★**** 

ADD TO QUEUE 

it-k'k'k’k'k'k'ki^'k-k-k'ki^'kif-k'k’kifikif’kitik'k'ki^'k’k/ 

/* This functions adds a new • node to the 
priority queue. It puts it in the last 
Dosition and then 'sifts it up' until it is 
in its proper location */ 
addq (ptr) struct trnode *ptr? < 

struct trnode *temp; 
i n t i f j ; 

/* Put the new node at the bottom */ 

qsize += 1? 

if (qsize == NUMQ) printf(fd,"Q OVERFLOWO); 
q(qsireJ = ptr; 

/* j is the node being sifted up */ 

i = qsize; 

/* i is the parent of j */ 

i = J/2; 

/* K'/hile j is not the root */ 

while(j > 1) { 

/* compare priorityofjtoi*/ 
i f ( q C i ] ->c os t < q(jJ->cost) return; 
else if (qlil->cost == qCjl->cost) ( 

if (q[i]->deoth >= q(j]->depth) return; 

> 

/* Exchange nodes and continue */ 

temp = q ( j 1 ; 
q t j 1 = q ( i 1 ; 
q ( i J = temp ; 
j = i ; 
i = j/2; 

> 

1 
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/★ ★ ★ ★ ★ 



GET TOP 

*****★************************/ 



/* This functions returns a pointer to the 
tree node which is on top of the priority 
queue. The queue is a heap. When the root 
of the full binary tree is removed^ the node 
in the last position is moved to the root and 
then 'sifted down' * / 
struct trnode *gettop() { 



struct trnode *saveot r , *t empJ 
i n t i / j ; 



savept r 

qtn = 

qsize = 

i = i; 

j = 2 * 

whi le( i 
i f 





/ * 


save 


off 


the 


root node */ 




= q 1 1 ) ; 


/* 


move 


the 


last 


node to the root 


★ / 


q (qs i ze) 


9 

/ 














/* 


dec remen t 


the 


qsize * / 




- i; 


/* 


sift 


the 


root 


node down */ 






/* 


i i s 


the 


parent node being looked 


i ; 


f* 


j i s 


the 


left 


ch i 1 d of i */ 




/* 


While 


were not out of the heap 


★/ 


<= qsize) { 











( j < qs i ze ) { 



/* Get the child of i with least cost */ 
if (qlj]->cost > q[j+l)->cost) j = j+l» 

/* if both the same cost/ get greater depth */ 
else if (q(jl->cost == q(j+ll->cost) 

i f ( q ( j } ->dePt h < q ( j + 1 1 ->aept h ) j = j + l? 

} 

/* If both children have lower priority/ quit */ 
if (qtiJ->cost < q(j)->cost) break; 
else if (q(i]->cost == q(jl->cost) 

if (q[il->depth >= qti]*>depth) break; 

/* exchange node i with j(its highest 
priority child) */ 

t emp = q C j ] ; 
q t j J = q [ i ) ; 
q [ i ) = temp ; 

/* set i = to j and do it again * / 

i = j / 
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4 



f| 




j 

> 



2 * i ? 



/* when done/ return the save off root */ 



ret urn ( savept r ) ; 
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APPENDIX B 



This appendix contains output listings from the various 
programs run to generate the plots described in Section IV. 
The first listing is the data for graph coloring without di- 
viding for graphs with from 10 to 30 vertices. The edge 
values are 30X» 50%, and 70%, 

The second listing is the output from the DIVIDE routine 
which divides the graphs and records data regarding the size 
of the biconnected components and the number of components. 
The graphs ranged in size from 6 to 100 vertices. 

The next listing contains data from tne routine to 
divide the graphs and color the components. Graphs with 25 
vertices and 2% to 1^% edges were used in this experiment. 

The last listing shows the effects of the amount of 
edges when n is held constant. For this experiment n = 20 
while the edges vary from 20% to 80% by 5% intervals. 
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COMPARISON OF 5 ORDERING METHODS FOR SEVERAL GRAPH SIZES 



A***************************** 



N = 10, 


EDGES = 


13(30%) 














ORDER 


BRANCHES 


NODES 
















avg 


max 


avg 


max 














RANDOM 


19 


65 


38 


1 15 














DEGREE 


10 


la 


21 


29 














OVEC 3 


9 


1 1 


21 


25 














ADJCEN 


9 


9 


21 


29 














LCA 


9 


9 


21 


29 


















BRANCH 


FACTOR BY 


LEVEL 










1 


2 


3 


9 


5 


6 


7 


8 


9 


10 


RANDOM 


1 .55 


1 .92 


1.96 


2.19 


1.75 


2.11 


1 .85 


2.93 


2.01 


0.00 


DEGREE 


1.23 


1.59 


1.82 


1 .99 


2.29 


2.91 


2.57 


2.78 


2.88 


0.00 


DVEC 3 


1 . 0 « 


1.95 


1.96 


2.22 


2.95 


2.97 


2.70 


2.83 


2.67 


0.00 


ADJCEN 


1.00 


1 .29 


2. 16 


2.99 


2.91 


2.95 


2,68 


2.63 


2.89 


0.00 


LCA 


1.00 


1.29 


2.15 


2.29 


2.97 


2.61 


2.79 


2.87 


2.92 


0.00 


N = 10, 


EDGES = 


22(50%) 














ORDER 


BRANCHES 


NODES 
















avg 


max 


avg 


max 














RANDOM 


29 


129 


53 


239 














DEGREE 


10 


18 


21 


35 














OVEC 3 


9 


1 1 


20 


25 














ADJCEN 


9 


1 1 


20 


25 














LCA 


9 


9 


20 


25 


















BRANCH 


FACTOR BY 


LEVEL 












1 


2 


3 


9 


5 


6 


7 


8 


9 


10 


RANDOM 


1.3a 


1.81 


1.72 


1 .99 


1.63 


2.29 


1 .57 


2.13 


1.75 


0.00 


DEGREE 


1.08 


1.36 


1.59 


1.76 


1 .93 


2.25 


2.57 


2.66 


3.13 


0.00 


OVEC 3 


1.01 


1.13 


1.92 


1.97 


2.27 


2.92 


2.56 


2.73 


3.13 


0.00 


ADJCEN 


1.00 


1.00 


1.35 


2.13 


2.99 


2.50 


2.65 


2.80 


3.16 


0.00 


LCA 


1.00 


1.00 


1.35 


2.07 


2.35 


2.95 


2.72 


2.87 


3.29 


0.00 


N = 10, 


EDGES = 


31(70%) 














ORDER 


BRANCHES 


NODES 
















avg 


max 


avg 


max 














RANDOM 


33 


160 


56 


250 














DEGREE 


10 


17 


19 


27 














OVEC 3 


9 


16 


18 


26 














ADJCEN 


9 


10 


18 


22 
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LCA 



9 



9 



18 



22 





1 


BRANCH FACTOR BY 
2 3 9 5 


LEVEL 

6 


7 


8 


9 


10 


RANDOM 


1.08 


1.50 


1.98 1.98 


1.51 


2.07 


1.93 


2.33 


1 .98 


0. 


DEGREE 


1.01 


1 .12 


1.29 1.36 


1.58 


1.98 


2.31 


2.63 


3.10 


0. 


OVEC 3 


1.02 


1.05 


1.19 1.21 


1.61 


2.15 


2.51 


2.73 


3.19 


0. 


ADJCEN 


1.00 


1.00 


1.00 1.13 


1.78 


2.91 


2.53 


2.72 


3.19 


0. 


LCA 


1.00 


1.00 


1.00 1 .13 


1.78 


2.23 


2.52 


2.76 


3.25 


0. 



N = 15, 


EDGES = 


31 (30%) 














ORDER 


BRANCHES 


NODES 
















avg 


max 


avg 


max 














RANDOM 


NOT 


ENOUGH SPACE AVAIL 












DEGREE 


22 


88 


98 


173 














DVEC 3 


17 


92 


38 


82 














ADJCEN 


15 


39 


36 


66 














LCA 


19 


22 


35 


99 


















BRANCH 


FACTOR BY 


LEVEL 












1 


2 


3 


9 


5 


6 


7 


8 


9 


10 


RANDOM 


NOT 


CALCULATED 














DEGREE 


1.21 


1 .59 


1.67 


1.86 


2.03 


2.03 


2.09 


2.22 


2.36 


2. 


DVEC 3 


1.06 


1.29 


1.77 


2.07 


2.11 


2.28 


2.31 


2.35 


2.53 


2. 


ADJCEN 


1.00 


1.02 


1.69 


2.13 


2.95 


2.98 


2.97 


2.96 


2.99 


2. 


LCA 


1.00 


1.02 


1.62 


2.07 


2.36 


2.91 


2.99 


2.57 


2.59 


2. 




1 1 


12 


13 


19 


15 












RANDOM 


NOT 


CALCULATED 














DEGREE 


2.66 


2.83 


3.05 


3.99 


0.00 


0.00 


0.00 


0.00 


0.00 


0. 


DVEC 3 


2.77 


2.88 


3.20 


3.93 


0.00 


0.00 


0.00 


0.00 


0.00 


0. 


ADJCEN 


2.78 


2.98 


3.15 


3.99 


0.00 


0.00 


0.00 


0.00 


0.00 


0. 


LCA 


2.82 


3.07 


3.32 


3.60 


0.00 


0.00 


0.00 


0.00 


0.00 


0. 



N = 15, 


EDGES = 


52(50%) 


ORDER 


BRANCHES 


NODES 




avg 


max 


avg 


max 


RANDOM 


NOT 


ENOUGH SPACE AVAIL 


DEGREE 


29 


72 


97 


129 


DVEC 3 


19 


97 


38 


89 


ADJCEN 


16 


61 


35 


115 


LCA 


15 


33 


32 


66 



loa 



00 

00 

00 

00 

00 

a9 

73 

65 

71 

00 

00 

00 

00 







BRANCH 


FACTOR BY 


LEVEL 










1 


2 


3 


a 


5 


6 


7 


8 


9 


10 


RANDOM 


NOT 


CALCULATED 














DEGREE 


l.ia 


1.30 


1.47 


1.57 


1.72 


1.81 


1.85 


2.07 


2.30 


2.25 


OVEC 3 


1 .04 


1.17 


1.34 


1 .59 


1.98 


1.9b 


2.04 


2.21 


2.22 


2.45 


ADJCEN 


1.00 


1.00 


1.04 


1.49 


2.02 


2.43 


2.39 


2.41 


2.25 


2.41 


LCA 


1.00 


1.00 


1.04 


1.47 


1.95 


2.21 


2.21 


2.29 


2.37 


2.50 




11 


12 


13 


14 


15 












RANDOM 


NOT 


CALCULATED 














DEGREE 


2.32 


2.57 


2.71 


3.14 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


OVEC 3 


2.49 


2.55 


2.70 


3.16 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


ADJCEN 


2.53 


2.63 


2.87 


3.22 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


LCA 


2.54 


2.70 


2.95 


3.38 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 



N = 15f 


EDGES = 


73(70X) 














ORDER 


BRANCHES 


NODES 
















avg 


max 


avg 


max 














RANDOM 


NOT 


ENOUGH SPACE AVAIL 












DEGREE 


24 


75 


44 


128 














OVEC 3 


21 


56 


38 


96 














ADJCEN 


15 


32 


31 


66 














LCA 


14 


28 


29 


55 


















BRANCH 


FACTOR BY 


LEVEL 












1 


2 


3 


a 


5 


6 


7 


8 


9 


10 


RANDOM 


NOT 


CALCULATED 














DEGREE 


1.02 


1.07 


1 .15 


1.31 


1.42 


1 .60 


1.71 


1.78 


1 .98 


2.20 


OVEC 3 


1.02 


1.02 


1.06 


1.24 


1.36 


1.60 


1.77 


2.01 


2.05 


2.33 


ADJCEN 


1.00 


1.00 


1.00 


1.00 


1.08 


1.49 


2.03 


2.31 


2.40 


2.49 


LCA 


1.00 


1.00 


1.00 


1.00 


1.08 


1.47 


1.93 


2.14 


2.14 


2.35 




1 1 


12 


13 


14 


15 












RANDOM 


NOT 


CALCULATED 














DEGREE 


2.21 


2.52 


2.70 


3.18 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


OVEC 3 


2.27 


2.46 


2.68 


3.15 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


ADJCEN 


2.45 


2.54 


2.85 


3.14 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


LCA 


2.49 


2.63 


2.91 


3.44 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 



N = 20 , EDGES = 57(30X) 



ORDER BRANCHES NODES 



avg max avg max 
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RANDOM 


NOT 


ENOUGH SPACE AVAIL 














DEGREE 


63 


312 


125 


601 
















OVEC 3 


33 


118 


70 


229 
















AOJCEN 


25 


208 


57 


418 
















LCA 


22 


121 


53 


247 




















BRANCH 


FACTOR BY 


LEVEL 












1 


2 


3 


4 


5 


6 


7 


8 


9 


10 




RANDOM 


NOT 


CALCULATED 
















DEGREE 


1.53 


1 .59 


1.72 


1.82 


1.86 


1.60 


1.76 


2.08 


2.08 


2. 


1 1 


DVEC 3 


1 .25 


1.34 


1 ,63 


1.87 


2.01 


2.05 


2.06 


2.28 


2.41 


2. 


32 


ADJCEN 


1.00 


1.00 


1.49 


1.84 


2,23 


2.22 


2.42 


2.40 


2.70 


2. 


76 


LCA 


1.00 


1.00 


1.49 


1.77 


2.07 


2.21 


2.36 


2.45 


2.51 


2. 


59 




1 1 


12 


13 


14 


15 


16 


17 


18 


19 


20 




RANDOM 


NOT 


CALCULATED 
















DEGREE 


2.17 


2.30 


2.41 


2.49 


2.51 


2.56 


2.76 


2.80 


3.39 


0. 


00 


DVEC 3 


2.36 


2.37 


2,38 


2.51 


2.57 


2.68 


2.79 


2.91 


3.26 


0. 


00 


ADJCEN 


2,70 


2.54 


2.38 


2.50 


2.55 


2.59 


2.81 


2,86 


3.31 


0. 


00 


LCA 


2.52 


2.54 


2.50 


2.56 


2.65 


2.78 


2.98 


3.18 


3.50 


0. 


00 



N = 20 , 


EDGES = 


95(50%) 
















ORDER 


BRANCHES 


NODES 


















avg 


max 


avg 


max 
















RANDOM 


NOT 


ENOUGH SPACE AVAIL 














DEGREE 


no 


1345 


203 


2592 
















DVEC 3 


76 


1206 


143 


2308 
















ADJCEN 


46 


881 


93 


1731 
















LCA 


31 


648 


65 


1281 




















BRANCH 


FACTOR BY 


LEVEL 














1 


2 


3 


4 


5 


6 


7 


8 


9 


10 




RANDOM 


NOT CALCULATED 
















DEGREE 


1 .15 


1.54 


1.57 


1.55 


1.79 


1 .93 


1.73 


1.59 


1 .64 


1 . 


71 


DVEC 3 


1.07 


1.26 


1.56 


1.60 


1.66 


1.70 


1.63 


1.59 


1.93 


1 . 


90 


AOJCEN 


1,00 


1.00 


1.00 


1.14 


1.68 


2.00 


2.1 1 


2.37 


2.25 


2. 


52 


LCA 


1.00 


1.00 


1.00 


1 .14 


1.65 


1.85 


1 .92 


2.15 


2.27 


2. 


32 




1 1 


12 


13 


14 


15 


16 


17 


18 


19 


20 




RANDOM 


NOT CALCULATED 
















DEGREE 


1 .89 


2.02 


2.01 


2.38 


2.41 


2.54 


2.66 


2.79 


3.09 


0. 


00 


DVEC 3 


2.13 


2.26 


2.15 


2.27 


2.45 


2.45 


2.66 


2.77 


3.05 


0. 


00 


ADJCEN 


2.39 


2.33 


2.30 


2.20 


2.40 


2.55 


2.61 


2.86 


3.1 1 


0. 


00 


LCA 


2.39 


2.36 


2.44 


2.33 


2.42 


2.55 


2.87 


3.09 


3.57 


0. 


00 
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I 

I 





N = 20, EDGES = I33(70r.) 



ORDER 


BRANCHES 


NODES 




avg 


max 


avg 


max 


RANDOM 


NOT 


ENOUGH 


SPACE AVAIL 


DEGREE 


52 


182 


91 


329 


DVEC 3 


36 


1 19 


67 


208 


ADJCEN 


21 




«3 


76 


LCA 


19 


35 


ai 


66 



BRANCH FACTOR BY LEVEL 





1 


2 


3 


a 


5 


6 


7 


8 


9 


10 


RANDOM 


NOT 


CALCULATED 














DEGREE 


1.02 


1.03 


1.07 


1.17 


1.45 


1.45 


1.67 


1.64 


1.88 


1.87 


DVEC 3 


1.00 


1.01 


1.05 


1.07 


1.25 


1 .56 


1.57 


1.75 


1 .82 


1 .52 


ADJCEN 1 


.00 1 


.00 1 


.00 1 


.00 1 


.00 1 


i.Ol 1 


.13 1 


i . 63 1 


1.97 2 


.30 


LCA 


1.00 


1 .00 


1 .00 


1.00 


1.00 


1 .01 


1 .12 


1.63 


1.88 


2 .oa 



RANDOM 


NOT 


CALCULATED 














DEGREE 


1.80 


1.67 


2.17 


2.29 


2.52 


2.64 


2.81 


3.10 


3.32 


0.00 


DVEC 3 


1.98 


2.11 


2.29 


2.43 


2.58 


2.72 


2.84 


3.07 


3.54 


0.00 


ADJCEN 


2.42 


2.52 


2.64 


2.80 


2.65 


2.76 


2.99 


3.13 


3.47 


0.00 


LCA 


2.23 


2.41 


2.53 


2.65 


2.77 


2.98 


3.19 


3.37 


3.91 


0.00 



N = 25, 


EDGES = 


90(30%) 
















ORDER 


BRANCHES 


NODES 


















avg 


max 


avg max 
















RANDOM 


NOT 


ENOUGH SPACE AVAIL 














DEGREE 


NOT 


ENOUGH SPACE AVAIL 














DVEC 3 


126 1158 


251 2288 
















ADJCEN 


63 


411 


132 819 
















LCA 


33 


115 


73 230 




















BRANCH FACTOR BY 


LEVEL 














1 


2 


3 4 


5 


6 


7 


8 


9 


10 




RANDOM 


NOT 


CALCULATED 
















DEGREE 


NOT 


CALCULATED 
















DVEC 3 


1.20 


1.35 


1.58 1.85 


1.80 


1 .93 


1 .98 


2.05 


1.79 


1 . 


84 


ADJCEN 


1.00 


1.00 


1.32 2.06 


2.12 


2.05 


2.20 


2.09 


2.15 


2. 


07 


LCA 


1.00 


1.00 


1.32 1.97 


1.91 


2.08 


1.98 


2.03 


2.08 


2. 


08 




1 1 


12 


13 14 


15 


lo 


17 


18 


19 


20 




RANDOM 


NOT 


CALCULATED 

















107 



degree not calculated 



OVEC 3 


1.77 


1.85 2.00 2.07 


2.05 


2.15 


2.32 


2.97 


2.37 


2.51 


AOJCEN 


2.09 


2.03 2.09 2.15 


2.26 


2.16 


2.32 


2.36 


2.52 


2.53 


LCA 


2.02 


2.03 2.09 2.22 


2.35 


2.92 


2.98 


2.59 


2.61 


2.75 




21 


22 23 2a 


25 












RANDOM 


NOT 


CALCULATED 














DEGREE 


NOT 


CALCULATED 














OVEC 3 


2.70 


2.78 2.85 3.97 


0.00 












AOJCEN 


2.79 


2.80 2.99 3.38 


0.00 












LCA 


2.89 


5.08 5.37 3.83 


0.00 













N 2 25, EDGES = 150(50%) 

ORDER BRANCHES NODES 

avg max avg max 

RANDOM NOT ENOUGH SPACE AVAIL 
DEGREE NOT ENOUGH SPACE AVAIL 
OVEC 3 178 1328 3a0 2568 

AOJCEN 97 932 193 1829 

LCA 38 139 76 271 

BRANCH FACTOR BY LEVEL 

123^56789 10 

RANDOM NOT CALCULATED 

DEGREE NOT CALCULATED 



OVEC 3 


1.09 


1 .29 


1.29 


1 .99 


1.76 


1.77 


1.73 


1 .83 


1.87 


1 .85 


AOJCEN 


1.00 


1.00 


1.00 


1.12 


1.90 


1 .93 


2.18 


2.26 


2.29 


2.1 1 


LCA 


1.00 


1 .00 


1.00 


1.12 


1.39 


1.86 


1.91 


2.08 


2.01 


1.99 





11 12 13 


19 


15 


16 


17 


18 


19 


20 




RANDOM 


NOT CALCULATED 
















DEGREE 


NOT CALCULATED 
















OVEC 3 


1.89 1.82 1.87 


1.80 


1.86 


1.92 


1.99 


2.03 


2.08 


2. 


16 


AOJCEN 


2.07 1.98 2.02 


1 .92 


2.09 


1 .97 


2.06 


2.01 


2. 12 


2. 


25 


LCA 


2.08 2.01 1.95 


1.95 


1.92 


2.09 


2.09 


2.20 


2.26 


2. 


37 





21 


22 23 


29 


25 




RANDOM 


NOT 


CALCULATED 






DEGREE 


NOT 


CALCULATED 






DVEC 3 


2.35 


2.59 2.66 


3.05 


0. 


00 


AOJCEN 


2.33 


2.38 2.59 


3.08 


0. 


00 


LCA 


2.59 


2.65 2.99 


3.29 


0. 


00 



108 



N = 25, EDGES = 210(70r.) 



ORDER 


BRANCHES NODES 








avg 


max avg max 






RANDOM 


NOT 


ENOUGH SPACE AVAIL 




DEGREE 


NOT 


ENOUGH SPACE AVAIL 




OVEC 3 


NOT 


ENOUGH SPACE AVAIL 




ADJCEN 


87 


723 169 1431 






LCA 


34 


122 64 225 










BRANCH FACTOR BY 


LEVEL 




1 


2 5 4 


5 


6 


RANDOM 


NOT 


CALCULATED 






DEGREE 


NOT 


CALCULATED 






DVEC 3 


NOT 


CALCULATED 






ADJCEN 


1.00 


1.00 1.00 1.00 


1.00 


1.01 


LCA 


1.00 


1.00 1.00 1.00 


1.00 


1.01 




11 


12 13 14 


15 


1 6 


RANDOM 


NOT 


CALCULATED 






DEGREE 


NOT 


CALCULATED 






OVEC 3 


NOT 


CALCULATED 






ADJCEN 


2.22 


2.34 2.22 2.16 


2.14 


2.13 


LCA 


1 .94 


2.05 2.12 1.95 


2.05 


2.02 




21 


22 23 24 


25 




RANDOM 


NOT 


CALCULATED 






DEGREE 


NOT 


CALCULATED 






OVEC 3 


NOT 


CALCULATED 






ADJCEN 


2.10 


2.44 2.59 2.85 


0.00 




LCA 


2.36 


2.58 2.84 3.41 


0.00 





7 8 9 10 



1.15 1 .as 1.90 2.06 
1.15 1 .a3 1.78 1.88 

17 18 19 20 



2. la 2.13 2.05 2.08 
2.08 2 . 0 a 2.10 2.27 



N = 30, EDGES = 130(30%) 



ORDER BRANCHES NODES 





avg 


max 


avg 


max 


RANDOM 


NOT 


ENOUGH 


SPACE A 


DEGREE 


198 


1504 


389 


2922 


OVEC 3 


88 


412 


179 


815 


ADJCEN 


63 


691 


134 


1391 


LCA 


31 


57 


71 


120 



BRANCH FACTOR BY LEVEL 
2 3 4 5 6 7 



1 



8 



9 



10 



RANDOM NOT CALCULATED 



DEGREE 


1.36 


1.73 


1 


.97 


1 


.02 


1.89 


1.97 


1.92 


1.86 


1.94 


1.93 


OVEC 3 


1.21 


1 .as 


1 


.64 


1 


.74 


1.94 


1.91 


1.92 


1.95 


1.92 


2.03 


ADJCEN 


1.00 


1.00 


1 


.30 


1 


.86 


2,16 


2.21 


2.26 


2.22 


2.32 


2.25 


LCA 


1.00 


1.00 


1 


.30 


1 


.83 


1 .94 


2.04 


2.10 


2.13 


2.13 


2. lo 





1 1 


12 


13 


14 


15 


lb 


17 


18 


19 


20 




RANDOM 


NOT 


CALCULATED 
















DEGREE 


2.09 


2.18 


2.00 


1.97 


2.05 


2.04 


2,12 


2.15 


2.37 


2. 


32 


OVEC 3 


1.95 


2.05 


2.16 


2.33 


2.25 


2.22 


2.24 


2.36 


2.30 


2. 


34 


ADJCEN 


2.28 


2.21 


2.25 


2.20 


2.31 


2.32 


2.21 


2.16 


2.31 


2. 


39 


LCA 


2.19 


2.14 


2.23 


2.19 


2.21 


2.21 


2.24 


2.32 


2.36 


2. 


38 





21 


22 


23 


24 


25 


26 


27 


28 


29 


30 


RANDOM 


NOT 


CALCULATED 














DEGREE 


2.26 


2.39 


2.43 


2.42 


2.51 


2.42 


2.59 


2.70 


2.80 


0.00 


OVEC 3 


2.41 


2.32 


2.54 


2.53 


2.62 


2.54 


2.59 


2.63 


2.84 


0.00 


ADJCEN 


2.41 


2.42 


2.42 


2.58 


2.51 


2.50 


2.55 


2.79 


2.87 


0.00 


LCA 


2.46 


2.51 


2.52 


2.60 


2.72 


2.81 


3.05 


3.12 


3.41 


0.00 



N = 30, 


EDGES = 217(50%) 














ORDER 


BRANCHES NODES 
















avg 


max avg max 




- 










RANDOM 


NOT 


ENOUGH SPACE AVAIL 












DEGREE 


NOT 


ENOUGH SPACE AVAIL 












OVEC 3 


289 


1825 555 3652 














ADJCEN 


79 


1995 163 3959 














LCA 


32 


76 70 148 


















BRANCH FACTOR 


BY 


LEVEL 












1 


2 3 4 


5 


6 


7 


8 


9 


10 


RANDOM 


NOT CALCULATED 














DEGREE 


NOT CALCULATED 














OVEC 3 


1.08 


1.12 1.33 1.59 1 


.70 


1.91 


1.89 


2.13 


1.99 


1.82 


ADJCEN 


1.00 


1.00 1.00 1.01 1 


.13 


1 .64 


2.01 


2.20 


2.25 


2.33 


LCA 


1.00 


1.00 1.00 1.01 1 


.12 


1.63 


1.88 


2.00 


2.01 


2.08 





11 12 13 14 


15 


16 


17 


18 


19 


20 


RANDOM 


NOT CALCULATED 














DEGREE 


NOT CALCULATED 














OVEC 3 


1.82 1,80 1,76 1.94 


1.89 


2.04 


2.12 


2.18 


2.17 


2.23 


ADJCEN 


2.30 2.29 2.33 2.27 


2.28 


2.32 


2.34 


2.13 


2.28 


2.39 



LCA 2.18 2.09 2.13 2.20 2.15 2,18 2.29 2,29 2.39 2.90 



21 22 23 29 

RANDOM NOT CALCULATED 

DEGREE NOT CALCULATED 

OVEC 3 2.33 2.93 2.99 2.51 

ADJCEN 2.90 2.37 2.93 2.95 

LCA 2.50 2.55 2.66 2.69 



25 


26 


27 


28 


29 


30 


2.59 


2.62 


2.83 


2.99 


3.17 


0.00 


2.58 


2.62 


2.78 


2.96 


3.25 


0.00 


2.80 


2.92 


3.06 


3.31 


3.75 


0.00 



N = 30, EDGES = 309(70X) 

ORDER BRANCHES NODES 

avg max avg max 



RANDOM NOT ENOUGH SPACE AVAIL 
DEGREE NOT ENOUGH SPACE AVAIL 
OVEC 3 215 1519 399 2896 

ADJCEN 96 1003 96 2003 

LCA 30 67 65 132 

BRANCH FACTOR BY LEVEL 



12 3 9 

RANDOM NOT CALCULATED 

DEGREE NOT CALCULATED 

OVEC 3 1.00 1,00 1.09 1.05 

ADJCEN 1.00 1.00 1.00 1.00 

LCA 1.00 1.00 1.00 1,00 

11 12 13 19 

RANDOM NOT CALCULATED 

DEGREE NOT CALCULATED 

DVEC 3 1.97 2.03 2.09 1.99 

ADJCEN 1.67 1.97 2.19 2.35 

LCA 1.66 1.88 1.96 2.09 

21 22 23 29 

RANDOM NOT CALCULATED 

DEGREE NOT CALCULATED 

DVEC 3 2.33 2.99 2.69 2.70 

ADJCEN 2.60 2.69 2.89 2.81 

LCA 2.68 2.70 2.86 2.98 



5 6 7 8 9 10 



1.08 1.22 1.36 1.39 1.57 1.65 

1.00 1.00 1.00 1.00 1.05 1.28 

1.00 1.00 1,00 1.00 1.05 1.28 

15 16 17 18 19 20 



1,83 1.72 1.68 2.07 2.16 2.29 
2.98 2.52 2.52 2.95 2.58 2.55 
2.10 2.26 2.39 2.59 2.95 2.95 

25 26 27 28 29 30 



2.89 2.98 3.15 3.39 3.75 0.00 
2.93 3.09 3.18 3.39 3.83 0.00 
5.13 3.28 3.98 3.75 9.23 0.00 



OUTPUT FROM GRAPH DIVIDING ROUTINE 
************************* 



1 1 1 



* * * 



N 



6 



8 



10 



12 



14 



EDGES AVG # SUBS AVG SIZE LARGEST SU 0 







n 


edges 


6 


3.4 


3.6 


3.6 


7 


2,5 


4.4 


5.2 


8 


1.7 


5.3 


7.1 


9 


1,3 


5.7 


8.7 


10 


1 . 1 


5.9 


9.9 


1 1 


1.0 


6.0 


11. 0 


12 


1.0 


6.0 


12.0 


13 


1.0 


6.0 


13.0 


0 


5.4 


3.5 


3.5 


10 


2.6 


6.1 


7.6 


12 


1.4 


7.5 


11.3 


14 


1 . 1 


7.9 


13.9 


16 


1,0 


8.0 


15.9 


18 


1.0 


8.0 


18.0 


20 


1.0 


0.0 


20.0 


22 


1.0 


8.0 


22.0 


10 


7,0 


4.0 


4.0 


12 


4.2 


6.3 


7.7 


14 


2.7 


0.2 


12.0 


16 


1.9 


9.0 


14.0 


18 


1.4 


9.5 


17.5 


20 


1.2 


9.7 


19.5 


22 


1 , 1 


9.9 


21.8 


24 


1.0 


10.0 


24.0 


26 


1.0 


10.0 


26.0 


28 


1.0 


10.0 


20.0 


12 


8.8 


4.2 


4.2 


15 


4.3 


0.1 


10.4 


10 


2.5 


10.2 


15.5 


21 


1.7 


11.3 


20.0 


24 


1.4 


1 1.6 


23.3 


27 


1 . 1 


11.9 


26.0 


30 


1.0 


12.0 


30.0 


33 


1.0 


12.0 


33.0 


36 


1.0 


12.0 


36.0 


39 


1.0 


12.0 


39.0 


14 


10.4 


4.6 


4.6 


17 


6.0 


0.1 


10.3 


20 


3.5 


11.1 


16.4 


23 


2.4 


12.5 


21 . 1 


26 


1 .8 


13.1 


24.8 



112 



29 


1.5 


13.5 


28.2 


32 


1.2 


13.8 


31 .5 


35 


1.2 


13.8 


34.6 


38 


1 . 1 


13.9 


37.9 


41 


1.0 


14.0 


41.0 


44 


1.0 


14.0 


44.0 


47 


1.0 


14.0 


47.0 


16 


12.9 


4.1 


4.1 


20 


5.3 


10.7 


13.8 


24 


1.5 


15.2 


22.8 


28 


1 . 1 


15.8 


27.7 


32 


1 . 1 


15.9 


31.9 


36 


1.0 


15.9 


36.0 


40 


1.0 


16.0 


40.0 


44 


1.0 


16.0 


44.0 


18 


14.1 


4.9 


4.9 


22 


7.4 


10.5 


13.4 


26 


4.4 


13.9 


20.9 


30 


2.8 


15.9 


27.2 


34 


2.1 


16.9 


32.7 


38 


1 .6 


17.2 


36.7 


42 


1.4 


17.6 


41.3 


46 


1.2 


17.7 


45.5 


50 


1.1 


17.9 


49.9 


54 


1.0 


17.9 


53.8 


58 


1.0 


18.0 


58.0 


62 


1.0 


18.0 


62.0 


66 


1.0 


18.0 


66.0 


70 


1.0 


18.0 


t > 9.9 


20 


16.4 


4.6 


4.6 


25 


7.5 


12.3 


16.0 


30 


4.1 


16.2 


24.7 


35 


2.4 


18.5 


32.6 


40 


1.8 


19.1 


38.8 


45 


1.4 


19.6 


44.3 


50 


1.2 


19.8 


49.8 


55 


1 . 1 


19.8 


54.8 


60 


1.0 


20.0 


59.8 


65 


1.0 


20.0 


65.0 


70 


1.0 


20.0 


69.9 


75 


1.0 


20.0 


75.0 


80 


1.0 


20.0 


80.0 


85 


1.0 


20.0 


85.0 


22 


18.0 


5.0 


5.0 



113 



27 


0.9 


12.9 


16.0 


32 


5.5 


16.5 


24.7 


37 


3.5 


18.9 


32.6 


a 2 


2.3 


20.5 


40.0 


a 7 


1.7 


21.1 


45.5 


52 


1 .5 


21.4 


50.6 


57 


1.3 


21.7 


56.1 


62 


1.2 


21.8 


61.5 


67 


1 .2 


21.0 


66.4 


72 


1 . 1 


21 .9 


71.8 


77 


1.0 


22.0 


77.0 


02 


1.0 


22.0 


02.0 


87 


1.0 


22.0 


07.0 


92 


1.0 


22.0 


92.0 


2 a 


19.0 


5.2 


5.2 


30 


8.4 


15.2 


19.7 


36 


4.5 


19.7 


30.1 


42 


2.6 


22.1 


39.3 


48 


1.7 


23.1 


46.4 


54 


1.3 


23.7 


53.4 


60 


1 .2 


23.8 


59.0 


66 


1 . 1 


23.9 


05.8 


72 


1.0 


24.0 


72.0 


78 


1.0 


24.0 


70.0 


84 


1.0 


24.0 


84.0 


PO 


1.0 


24.0 


90.0 


26 


21.0 


5.2 


5.2 


32 


10.5 


15.1 


19.8 


38 


6 . 1 


19.9 


30.0 


44 


3.9 


22.3 


30.5 


50 


2.7 


23.8 


46.5 


56 


2 . 1 


24.7 


53.6 


62 


1.6 


25.2 


60.6 


60 


1.3 


25.6 


67.1 


74 


1.2 


25.7 


73.0 


00 


1.1 


25.9 


79.7 


06 


1.0 


26.0 


85.9 


92 


1.0 


25.9 


91 .6 


96 


1.0 


26.0 


98.0 


104 


1.0 


26.0 


104.0 


no 


1.0 


26.0 


110.0 


116 


1.0 


26.0 


116.0 


28 


23.6 


5.4 


5,4 


35 


10.5 


17.0 


22.4 


42 


5.5 


22.8 


35.5 



1 la 



49 


3,4 


25.3 


45.5 


56 


2.2 


26.6 


54.1 


65 


1,6 


27.2 


01.3 


70 


1.3 


27.6 


69.1 


77 


1.2 


27.7 


76.3 


84 


1.1 


27.9 


83.8 


91 


1 . 1 


27.9 


90.9 


98 


1.0 


28,0 


98.0 


105 


1.0 


28.0 


104.8 


112 


1,0 


28.0 


112.0 


119 


1.0 


28.0 


1 19.0 


30 


25.5 


5,6 


5.6 


37 


1 1.9 


17.6 


22.7 


44 


6.7 


23.5 


36.0 


51 


4.4 


25.7 


44.8 


58 


2.7 


28.1 


55.6 


65 


2.0 


28.7 


63.0 


72 


1.5 


29.3 


70.3 


79 


1.5 


29.5 


78.2 


86 


1.3 


29.7 


85.5 


93 


1.1 


29.9 


92.8 


100 


1,0 


29.9 


99.7 


107 


1.0 


30.0 


107.0 


114 


1.0 


30.0 


114.0 


121 


1.0 


30.0 


121.0 


128 


1.0 


30.0 


128.0 


135 


l.O 


30.0 


135.0 


40 


34,8 


6.2 


6.2 


60 


7.9 


31.7 


49.0 


80 


2.7 


37,8 


75.9 


100 


1.5 


39.4 


98.6 


120 


1 . 1 


39,9 


119.8 


140 


1.0 


40.0 


140.0 


160 


1.0 


40,0 


160.0 


180 


1.0 


40.0 


180.0 


50 


44.5 


6.5 


6.5 


75 


10.6 


38.7 


61.2 


100 


4.0 


46.7 


95.6 


125 


2.0 


48.5 


120.8 


150 


1.4 


49.5 


148.7 


175 


1 . 1 


49,9 


1 74.9 


200 


1.0 


50.0 


200,0 


225 


1.0 


50.0 


225,0 


250 


1.0 


50.0 


250.0 



115 



60 


54.3 


6.7 


6.7 


90 


12.7 


46.1 


72.2 


120 


4.5 


56.0 


113.9 


150 


2.1 


58.7 


147.5 


180 


1.2 


59.7 


179.0 


210 


1.1 


59.9 


209.6 


2ao 


1.0 


60.0 


240.0 


270 


1.0 


60.0 


270.0 


iOO 


1.0 


60.0 


300.0 


530 


1.0 


60.0 


330.0 


70 


64.3 


6.7 


6.7 


105 


15.1 


53.7 


85.0 


lao 


5.7 


64.6 


132.1 


175 


2.6 


68.2 


171.9 


210 


1.5 


69.3 


207.8 


295 


1.2 


69.8 


244.8 


280 


1.0 


70.0 


280.0 


315 


1.0 


70.0 


315.0 


350 


1.0 


70.0 


350.0 


80 


74.3 


6.7 


6.7 


120 


14.6 


63.9 


99.8 


160 


4.7 


75.8 


153.6 


200 


2.0 


78.7 


197.7 


240 


1.3 


79.7 


239.6 


280 


1.0 


80.0 


280.0 


320 


1.0 


80.0 


320.0 


360 


1.0 


80.0 


360.0 


400 


1.0 


80.0 


400.0 


90 


83.8 


7.2 


7.2 


155 


19.1 


68.6 


107.8 


180 


6.5 


83.4 


lo 9.9 


225 


3.1 


87.6 


220.8 


270 


1.9 


89.0 


267.4 


315 


1.2 


89.8 


314.5 


360 


1.0 


89.9 


3 o 0.0 


405 


1.0 


90.0 


405.0 


450 


1.0 


90.0 


450.0 


495 


1.0 


90.0 


495.0 


100 


93.0 


8.0 


8.0 


150 


20.9 


77.1 


122.5 


200 


7.2 


91.9 


186.0 


250 


3.0 


97.8 


246.4 


500 


1.6 


99.4 


299.1 


350 


1.2 


99.7 


349.0 



116 



★ * 



400 


1 . 1 


99.9 


399.2 


450 


1 .0 


100.0 


450.0 


500 


1.0 


100.0 


500.0 


550 


1.0, 


100.0 


550.0 



OUTPUT FROM DIVIDE AND COLOR PROGRAM 



N = as, EDGES = 6( 2X) 



ORDER 


BRANCHES 


NODES 




avg max 


avg 


max 


RANDOM 


6 


6 


12 


12 


DEGREE 


6 


6 


12 


12 


DEG 3 


6 


6 


12 


12 


ADJCEN 


6 


6 


12 


12 


LCA 


6 


6 


12 


12 


N = 25, 


EDGES 


- 


12( 


4%) 



ORDER 


BRANCHES 


NODES 




avg max 


avg 


max 


RANDOM 


12 


12 


23 


24 


DEGREE 


12 


12 


23 


24 


DEG 3 


12 


12 


23 


24 


ADJCEN 


12 


12 


23 


24 


LCA 


12 


12 


23 


24 


N = 25, 


EDGES 




10( ( 


b%) 



ORDER 


BRANCHES 


NODES 




avg 


max 


avg 


max 


RANDOM 


17 


18 


34 


36 


DEGREE 


17 


19 


34 


36 


DEG 3 


17 


18 


34 


36 


ADJCEN 


17 


18 


34 


36 


LCA 


17 


18 


34 


36 



N = 25, 


EDGES = 


24( 8X) 


ORDER 


BRANCHES 


NODES 




avg 


max 


avg max 


RANDOM 


21 


25 


43 55 



* 



117 



I 




DEGREE 


22 


a6 


46 


90 


DEG 3 


21 


27 


44 


55 


ADJCEN 


21 


23 


45 


54 


LCA 


21 


23 


44 


54 



N = 25, 


EDGES = 


30(10%) 


ORDER 


BRANCHES 


NODES 




avg 


max 


avg 


max 


RANDOM 


25 


45 


53 


91 


DEGREE 


27 


52 


57 


107 


DEG 3 


24 


31 


52 


68 


ADJCEN 


25 


30 


52 


67 


LCA 


23 


30 


52 


68 



N = 25, 


EDGES = 


36(12%) 


ORDER 


BRANCHES 


NODES 




avg 


max 


avg 


max 


RANDOM 


51 


945 


108 


1893 


DEGREE 


34 


191 


72 


380 


DEG 3 


26 


94 


59 


189 


ADJCEN 


24 


35 


57 


82 


LCA 


24 


35 


57 


82 



N = 25, 


EDGES = 


42(14%) 


ORDER 


BRANCHES 


NODES 




avg 


max 


avg 


max 


RANDOM 


84 


1012 


172 


2021 


DEGREE 


43 


208 


91 


425 


DEG 3 


29 


86 


64 


177 


ADJCEN 


29 


170 


65 


349 


LCA 


25 


54 


58 


109 



* 



* 



************************ 
RESULTS OF VARYING EDGES FOR A CONSTANT N 



N = 20, EDGES = 58(20%) 



order branches 
avg max 



nodes 
avg max 



Loga r ithm of max 
branches nodes 
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DEGREE 


46 


363 


96 


731 


5.89 


6.59 


DVEC 3 


27 


137 


60 


278 


4.92 


5.63 


ADJCEN 


23 


82 


52 


166 


4.41 


5.1 1 


LCA 


20 


33 


46 


70 


3.50 


4.25 



N = 20 , 


EDGES = 


47(25%) 






order 


branches 


nodes 


Logar i t hm 


of max 




avg 


max 


avg 


max 


branches 


nodes 


DEGREE 


45 


127 


94 


248 


4.84 


5.51 


DVEC 3 


28 


93 


62 


182 


4.53 


5.20 


ADJCEN 


21 


43 


52 


87 


3.76 


4.47 


LCA 


20 


27 


48 


63 


3.30 


4.14 



N = 20, 


EDGES = 


57(30%) 






order 


branches 


nodes 


Logar i t hm 


of max 




avg 


max 


avg 


max 


branches 


nodes 


DEGREE 


49 


187 


99 


380 


5.23 


5.94 


DVEC 3 


32 


164 


67 


320 


5.10 


5.77 


ADJCEN 


25 


109 


56 


223 


4.69 


5.41 


LCA 


21 


55 


48 


117 


4.01 


4.76 



N = 20, 


EDGES = 


66(35%) 






order 


branches 


nodes 


Logar i t hm 


of max 




avg 


max 


avg 


max 


branches 


nodes 


DEGREE 


112 


1 119 


219 


2198 


7.02 


7.70 


DVEC 3 


52 


308 


104 


608 


5.73 


6.41 


ADJCEN 


38 


513 


80 


1015 


6.24 


6.92 


LCA 


23 


62 


51 


123 


4.13 


4.81 



N = 20, 


EDGES = 


76(40%) 






order 


branches 


nodes 


Logar i t hm 


of max 




avg 


max 


avg 


max 


branches 


nodes 


DEGREE 


NOT 


ENOUGH SPACE AVAIL 




DVEC 3 


59 


267 


116 


524 


5.59 


6.26 


ADJCEN 


30 


292 


65 


593 


5.68 


6.39 


LCA 


23 


285 


52 


568 


5.65 


6.34 


N = 20, 


EDGES = 


85(45%) 
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order 


branches 


nodes 


Loga r i t hm 


of max 




a vg 


max 


a vq 


max 


branches 


nodes 


DEGREE 


92 


635 


173 


1206 


6.45 


7.10 


OVEC 3 


56 


425 


109 


818 


6.05 


6.71 


ADJCEN 


37 


379 


77 


749 


5.94 


6.62 


LCA 


23 


64 


49 


128 


4.16 


4.85 



Z 

II 

ru 

o 


EDGES = 


95(50%) 






order 


branches 


nodes 


Logar 1 1 hm 


of max 




avg 


max 


avg 


max 


branches 


nodes 


DEGREE 


101 


811 


188 


1608 


6.70 


7.38 


DVEC 3 


57 


367 


107 


676 


5.91 


6.52 


ADJCEN 


34 


501 


71 


980 


6.22 


6.89 


LCA 


28 


584 


60 


1150 


6.37 


7.05 



N = 


EDGES = 


104(55%) 






order 


branches 


nodes 


Loga r i t hm 


of max 




avg 


max 


avg 


max 


branches 


nodes 


DEGREE 


106 


431 


194 


777 


6.07 


6.66 


DVEC 3 


69 


440 


131 


844 


6.09 


6.74 


ADJCEN 


34 


258 


70 


497 


5.55 


6.21 


LCA 


23 


79 


48 


156 


4.37 


5.05 



N = 20 , 


EDGES = 


114(60%) 






order 


branches 


nodes 


Loga r i t hm 


of max 




avg 


max 


avg 


max 


branches 


nodes 


DEGREE 


99 


9b9 


130 


1696 


6.38 


7.44 


OVEC 3 


62 


576 


114 


1029 


6.36 


6.94 


ADJCEN 


37 


561 


75 


1088 


6.33 


6.99 


LCA 


27 


438 


56 


855 


6.08 


6.75 


N = 20 , 


EDGES = 


123(65%) 






order 


branches 


nodes 


Loga r i t hm 


of max 




avg 


max 


avg 


max 


branches 


nodes 


DEGREE 


86 


765 


155 


1351 


6.64 


7.21 


DVEC 3 


55 


468 


101 


832 


6.15 


6.72 


ADJCEN 


28 


188 


57 


365 


5.24 


5.90 


LCA 


21 


53 


43 


103 


3.97 


4.63 


N = 20, 


EDGES = 


133(70%) 
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orde r 


branches 


nodes 


Logari thm 


of max 




avg 


max 


avg 


max 


branches 


nodes 


DEGREE 




419 


130 


748 


6.04 


6.62 


DVEC 5 


47 


282 


86 


498 


5.64 


6.21 


ADJCEN 


21 


57 


44 


114 


4.04 


4.74 


LCA 


20 


33 


41 


62 


3.50 


4.13 



H = 20, EDGES = 1^2(75X3 



order 


branches 


nodes 


Logari thm 


of max 




avg 


max 


avg 


max 


branches 


nodes 


DEGREE 


60 


295 


103 


499 


5.69 


6.21 


DVEC 3 


38 


226 


69 


403 


5.42 


6.00 


ADJCEN 


22 


66 


43 


118 


4.19 


4.77 


LCA 


19 


37 


39 


71 


3.61 


4.26 



N = 20, EDGES = 152(80X) 



order 


branches 


nodes 


Logari thm 


of max 




avg 


max 


avg 


max 


branches 


nodes 


DEGREE 


45 


143 


76 


237 


4.96 


5.47 


DVEC 3 


30 


71 


54 


118 


4.26 


4.77 


ADJCEN 


20 


46 


39 


81 


3.83 


4.39 


LCA 


19 


27 


37 


49 


3.30 


3.89 
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